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Abstract 



The Higgs low-energy theorem gives a simple and elegant way to estimate the couplings 
of the Higgs boson to massless gluons and photons induced by loops of heavy parti- 
cles. We extend this theorem to take into account possible nonlinear Higgs interactions 
resulting from a strong dynamics at the origin of the breaking of the electroweak sym- 
• • metry. We show that, while it approximates with an accuracy of order a few percents 

. £h single Higgs production, it receives corrections of order 50% for double Higgs produc- 

tion. A full one-loop computation of the gg — > hh cross section is explicitly performed 
in MCHM5, the minimal composite Higgs model based on the SO(5)/SO(4) coset with 
the Standard Model fermions embedded into the fundamental representation of SO (5). 
In particular we take into account the contributions of all fermionic resonances, which 
give sizeable (negative) corrections to the result obtained considering only the Higgs 
nonlinearities. Constraints from electroweak precision and flavor data on the top part- 
ners are analyzed in detail, as well as direct searches at the LHC for these new fermions 
called to play a crucial role in the electroweak symmetry breaking dynamics. 



1 Introduction 



There is growing evidence that a Higgs boson is the agent of electroweak symmetry breaking 
(EWSB) |lj. And one of the most pressing questions is to uncover the true nature of this Higgs 
boson: Is it the elementary scalar field of the Standard Model? Is it part of a supermultiplet? 
Is it a composite scalar emerging as a bound state from a strongly coupled sector? Theoretical 
arguments based on naturalness considerations tend to favor one of the two later scenarios. But the 
LHC experiments have now opened a new data-driven era and the first experimental indications 
might come from possible deviations in the measurements of the Higgs couplings [2] compared to 
the ones predicted by the Standard Model (SM) that are unambiguously fixed by the value of 
the Higgs mass itself. But both supersymmetric and composite Higgs bosons in the decoupling 
limits can be arbitrarily close to the SM Higgs at the energy scale currently probed by the LHC 
and deciphering the different scenarios from one another might require more luminosity than the 
one currently accumulated. A more unambiguous answer would come from a direct observation 
of additional particles: supersymmetric partners of the SM particles or additional resonances of 
the strong sector. In both scenarios, new particles in the top sector have a special status: in 
supersymmetric models, the stops cannot be too heavy without destabilizing the weak scale (31; the 
top partners in the composite models are responsible for generating the potential of the would-be 
Goldstone Higgs boson and, as recently noticed in Refs. (4j, they have to be lighter than about 
700 GeV to accommodate a Higgs boson as light as 125 GeV. Even if the actual numbers are 
model-dependent, these conclusions are rather generic and they certainly call for improving the 
ongoing dedicated direct searches for the top partners (H]. 

In principle, indirect information on these top partners could also be obtained in Higgs physics. 
It is indeed quite ironic that, while the Higgs boson is supposed to be at the origin of the masses 
of all elementary particles, the current most sensitive channels are the ones that involve massless 
particles, i.e. particles with no direct coupling to the Higgs boson: the gluons for its production 
and the photons for its decay. Clearly within the SM, these processes appear at the loop level 
only and are therefore potentially sensitive to new states circulating in the loops. The structure 
of the Higgs couplings to photons and gluons are beautifully captured by the background field 
method known as Low-Energy Theorem (LET) [6j[7] and the effects of the top partners, or any 
new particles, are encoded by their sole contributions to the QED and QCD beta functions. In 
composite Higgs models, the corrections due to the Higgs compositeness to the SM gluon fusion 
production cross section, a(gg — > h), and to the SM decay width into two photons, F(h — > 77), 
were estimated in Ref. [8] to be generically of the order of £ = (v/f) 2 , where v ~ 246 GeV is the 
weak scale and / is the characteristic scale of the strong sector (the equivalent of f n in QCD). From 
dimensional arguments, the additional corrections due to the top partners should be of the order of 
(mt/rriT) 2 , TTit T being respectively the mass of the top and the typical mass scale of his partners, 
i.e. a correction potentially as big as the one originating from the strong dynamics itself when the 
top partners are around 700 GeV. However, it was quickly realized [9-11 that these leading-order 
corrections from the top partners actually cancel and that they only give a contribution that we 
shall estimate (see Section 3.3 ) to scale like ^{frih/mt) 2 /10, i.e. one order of magnitude smaller than 
the strong dynamics contribution, sweeping any hope to learn anything about the top partners from 
the measurement of the Higgs production by gluon fusion and leaving only the Higgs production 
in association with a top-antitop pair [12] as a place to indirectly look for new physics in the first 
LHC run. A proper effective description of the top partners/fermionic resonances would be needed 
to study this promising channel. 

The second LHC run will certainly increase the sensitivity of the direct searches for top partners. 
But it will also open new possibilities to indirectly probe the top partner sector in Higgs physics by 
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exploring multi-Higgs production. A generalization of the LET exists for double-Higgs production 
in the SM, and it is known to give a reasonably good estimate of the total rate within 20% 



accuracy for a light Higgs 13,14 but it badly fails to reproduce the differential distributions 15 



The extension of the LET to strong EWSB models is not totally trivial since, in the SM, there is a 
strong destructive interference between two contributions and therefore strong dynamics that gives 
rise to a third new contribution can have order-one effects on the double Higgs production by gluon 



fusion, gg — > hh, as already noticed in Refs. 16-18 . In an explicit composite Higgs model, we are 



going to compare the LET results to an explicit one-loop computation taking into account all the 
contributions from the fermionic resonances in the top sector. On the other hand the pure strong 
dynamics effects arising due to the sole Higgs nonlinearities are overestimating the production rate 
which receives 0(30%) corrections from the top partners. 

The outline of our paper is as follows. In Section [2j we first present two equivalent effective 
Lagrangians, the linearly realized strongly-interacting light Higgs Lagrangian [8] as well as a chiral 
Lagrangian where the SU(2)l x U{\)y SM gauge symmetry is nonlinearly realized 19 , to describe 
Higgs physics at the LHC and we relate these effective Lagrangians to explicit composite Higgs 
models. In Section |3j we recall the LET in the SM and extend it to composite Higgs models, repro- 
ducing the often heard results that the gluon fusion production cross sections can be significantly 
reduced by up to 50%-^ 70% compared to the SM production. We also re-discuss within the LET 
approximation the cancellation of the top partner contributions to single Higgs production by gluon 
fusion and we estimate the size of the corrections to the LET results. We finally derive the LET 
for double Higgs production by gluon fusion and show that a factor 2-^3 enhancement over the SM 
can be obtained in specific examples. This LET enhancement factor is nonetheless smaller than 
the one obtained considering only the pure Higgs nonlinearities effects. Section [4] is devoted to the 
top partners in an explicit minimal composite Higgs model, MCHM5, based on the coset space 
SO(5)/ SO(4) with SM fermions embedded into fundamental representations of SO(5). We first 
discuss in detail the constraints on the masses and couplings of these new fermions coming from 
electroweak (EW) precision measurements and from flavor physics. We then present exhaustive 
constraints from direct searches for these top partners using the current Tevatron and LHC data. 
In Section [5j a detailed computation of the contribution of the top partners to a{gg — > h) is pre- 
sented, confirming that the LET gives a good approximation and confirming the scaling behavior 
of the corrections obtained in Section [3| Finally, Section [6] is devoted to double Higgs production. 
We demonstrate explicitly that the full loop computation with the top partners can easily give 
negative corrections of order 30% to the pure strong dynamics computation where only the Higgs 



nonlinearities were taken into account 17 . Contrary to the single Higgs production case, the LET 
does not capture well these contributions of the top partners, and this discrepancy is further am- 
plified when looking at particular regions of the phase space selected by kinematical cuts. Various 
technical details are collected in a series of Appendices. 



2 Low-energy effective Lagrangian for a composite Higgs boson 

An interesting solution to the hierarchy problem is given by the Higgs boson being a composite 
bound state emerging from a new strongly-interacting sector, broadly characterized by a mass 
scale m p and a coupling g p . If in addition the Higgs emerges as a pseudo-Goldstone boson of a 
spontaneous symmetry breaking Q /% at the scale / = m p /g p , then it can be naturally lighter than 
the other resonances of the strong sector and v 2 / f 2 <C 1 can be accommodated. A low-energy, 
model-independent description of this idea is given by the strongly-interacting light Higgs (SILH) 
Lagrangian (§], which applies to the general scenario where the Higgs is a light pseudo-Goldstone 
boson, including Little Higgs and Holographic composite Higgs models. At scales much smaller 
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than m p , deviations from the SM are parameterized in terms of a set of dimension-six operators. 
The subset that will be relevant in our discussion is 
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where g s , g, g' are the SM SU(3) C x SU(2) lxU (l)y gauge couplings, whereas A and are the Higgs 
quartic and Yukawa coupling appearing in the SM Lagrangian, respectively. The first four operators 
are genuinely sensitive to the strong interaction, whereas the last two parameterize the effective 
couplings of the Higgs to gluons and to photons, respectively, mediated by loops of heavy particles. 
As the operators proportional to c g and Cy do not respect the symmetry under which the Goldstone 
Higgs shifts they have to be suppressed by powers of the couplings which break this symmetry, 
thus explaining the extra factor g\^jg\ (gsM = Ut,g) appearing in front of these operators. These 
operators are important in the presence of relatively light resonances, i.e. when g p ~ gsM- Recent 
analyses |4| show that in a large class of models, a Higgs as light as 125 GeV implies the presence of 
one or more anomalously light (sub-TeV) fermionic resonances, which thus can contribute sizeably 
to c g and c 7 . Also notice that by choosing flavor- diagonal couplings yj Minimal Flavor Violation 
(MFV) is automatically implemented in the Lagrangian so that it complies with flavor bounds. 

In Eq. ([!]) we have kept explicitly the operator proportional to c r , which can be eliminated at 
0(1/ f 2 ) by a field redefinition 

H H + a(H^H)H/f , (2) 

under which 

ch — ^ ch + 2a , c r — > c r + 4a , cq — > cq + 4a , c y — )• c y — a , (3) 

while c g and c 7 do not change under this transformation. The choice c r = corresponds to the 
'SILH basis', which can be reached starting from a generic basis where c r ^ by applying the 
transformation in Eq. ^ with a = —c T j\. We choose to keep explicitly the operator proportional 
to c r as the 'natural' basis for nonlinear u-models actually corresponds to a non- vanishing c r |10|. 
Furthermore, since physical amplitudes have to be invariant under field redefinitions, Eqs. ^ and 
([3]) will be used as a consistency check of our results. In Eq. ([!]) we have also omitted the custodial 
breaking operator 

0^(VW) 2 , (4) 

where H^D^H = H^(D fl H) — (D IJ 'H)^H, which gives a contribution to the T parameter T = ctv 2 j f 2 
and thus is strongly constrained by electroweak data. This operator does not contribute to the pro- 
cesses that we will be interested in. If the strong sector is invariant under custodial symmetry, as it 
happens for example in models based on the coset SO(5) / SO(A), then ct vanishes at tree-level. For 
a discussion of Higgs physics where the assumption of custodial invariance is relaxed, see Ref. [20] . 

The SILH Lagrangian represents an expansion in £ = (v/ f) 2 and can be used in the vicinity of 
the SM, which corresponds to £ = 0. On the other hand, the technicolor limit (£ —> 1) requires the 
resummation of the full series in £. Such a resummation is possible in the Holographic Higgs models 



of Refs. [21-23 . These models are based on a five-dimensional gauge theory in Anti-de-Sitter (AdS) 



space-time. In the minimal realization, the bulk symmetry SO(5) x U(l) is broken to the SM group 
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SU(2) L x U(1) Y on the ultraviolet (UV) brane and to 50(4) x U(l) on the infrared (IR) brane. 
The coset SO(5)/SO(4) provides four Goldstone bosons, one of which is the physical Higgs boson 
and the three remaining ones are eaten by the massive SM vector bosons. The Higgs couplings 
to gauge bosons and its self-interactions are modified compared to the SM, and the modification 
factors can be expressed in terms of the parameter £. The Higgs Yukawa couplings and the form 
of the Higgs potential of the low-energy effective theory depend on the way the SM fermions are 
embedded into representations of the bulk symmetry. In the second part of this work we refer to 
the model MCHM5 [22] where the fermions transform in the fundamental representation of SO(5). 
An alternative realization of the SO(5)/SO(4:) composite Higgs, denoted by MCHM4, contains 
fermions embedded into the spinorial representation |21| (for more details see App.jB]). In this case, 
however, large corrections to the Zhj^bi coupling are present and rule out an important part of 
the parameter space ['2 lj. In contrast, if fermions are embedded into the fundamental or adjoint 
representation of SO (5), the custodial symmetry of the strong sector includes a left-right parity, 
which protects the Zbibi coupling from receiving tree- level corrections [25] . 

Another useful description of the low-energy theory is given by an effective chiral Lagrangian 
where the SU(2) x U(l)y symmetry is nonlinearly realized. The Goldstone bosons 7r a (a = 1,2,3) 
providing the longitudinal degrees of freedom of the and Z bosons are introduced by means of 
the field 



£(aO = e laana{x)/v 



(5) 



where v ~ 246 GeV and a a are the Pauli matrices. The field £ transforms linearly under SU(2)l x 
SU(2)r. Introducing a scalar field h, assumed to transform as a singlet under the custodial sym- 
metry, leads to the following effective Lagrangian 19 



l 2 (d»h?-V(h) + V -Tr 
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with the mass of the scalar given by m^. In Eq. ^ we have introduced the higher-dimensional 
couplings k g ,k2g,k~ / , which are mediated at loop level by strong sector resonances. The Higgs 
couplings to fermions, c, c%, are assumed to be flavor- diagonal, so that MFV is realized. In 
Table [l] the values of the couplings in the effective Lagrangian Eq. ^ are listed in the SILH 
approach and in the holographic Higgs model MCHM5. The SM with an elementary Higgs boson 



corresponds to a 
terms in h. 



ds = d 4 



1, c 2 



fog = k*y = and vanishing higher order 



3 Applying the Higgs Low-energy Theorem 

In this section we discuss applications of the Higgs Low-energy Theorem [6j[7] in composite models. 
The LET allows one to obtain the leading interactions of the Higgs boson with gluons and photons 
arising from loops of heavy particles. By heavy particles we mean here both SM states (W and top) 
and new states belonging to the composite sector. These couplings are needed in the computation 
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Parameters SILH MCHM5, no light resonances 
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Table 1: Values of the couplings of the effective Lagrangian Eq. (|6| in the SILH framework (with ct = 0) and 
for MCHM5 with no anomalously light fermionic resonances (the latter are taken from Ref. [26] )• The values 
of the SILH parameters in MCHM5 are, in the 'natural' basis for the nonlinear u-model where c r = —4 ch 
and neglecting effects of fermion resonances, ch = 1/3, c r = —4/3, c y = 4/3, c§ = —4/3 . 

of the cross sections of single and Higgs pair production via gluon fusion at the LHC as well as of 
the partial width of the decay h — > 77 . We will adopt a model-independent approach and compute 
these quantities in terms of the parameters of the effective Lagrangians defined in Section [2j Eqs. ([!]) 
and ([6]), putting special emphasis on the former, namely the SILH description. Our analysis extends 
the results of Refs. pTo|[2T] to Higgs pair production in gluon fusion, and also includes a discussion 
of corrections to the LET approximation arising from higher order terms in the 1/M expansion, 
where M is the mass of the generic heavy particle running in the loops. Notice that the LET can 
be extended to 2-loop order to include the leading QCD corrections, see for example Ref. [7j for 
applications in the SM. However, our discussion will be mainly limited to couplings at the leading 
1-loop order. 

3.1 Higgs interactions with gluons 

According to the LET the interactions of the physical Higgs boson with gluons, mediated by loops 
of heavy coloured particles, can be obtained by treating the Higgs H as a background field and 
taking the field-dependent mass of each heavy particle as a threshold for the running of the QCD 
gauge coupling]^] Assuming the heavy particles to transform in the fundamental representation of 
SU(3) C one obtains the following effective Lagrangian 

= S-^ Ga " V £ Sb * log < W . (7) 

Pi 

where 5b = 2/3 if particle pi is a Dirac fermion, and 5b = 1/6 if it is a complex scalar. In this paper 
we will focus only on the effects of the heavy fermion sector, which in composite Higgs models 

throughout the paper, we will denote by H both the Higgs doublet and the scalar field with {H} 7^ , as it will 
always be clear from the context which one we are referring to. On the other hand, h denotes the physical Higgs 
field. 
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typically includes new states beyond the top quark. By expanding the field-dependent masses of 
the heavy particles around the vacuum expectation value (VEV) (H) we obtain the couplings of 
the Higgs boson to gluons mediated by loops of heavy fermions 



99 



Q 2 

Us r~ia /-ia[iv 

96vr 2 



A x h + -A 2 h 2 + . . . 



where we have defined 



^logdet^ 2 (F; 



(H) 



(8) 



(9) 



with M 2 = Ai'M, and A4 is the heavy fermion mass matrix. In the SM only the top quark 
contribute^] with mt(H) = y t H/\/2, so that Eq.Q can be rewritten at all orders in h as (see for 
example Refs. [Tip]) 



(10) 



The corresponding gauge invariant operator is log(iJT H)G^ u G a ^ v , which is associated with a chiral 
fermion. The lowest-order operator arising from vector-like fermions is instead HG a ^ u G atlu . The 
effects of these two operators on double Higgs production were discussed in Ref. [29] . 

Using Eq. Q it is straightforward to derive the expression of the hgg and hhgg couplings in 
the SILH formalism. We refer the reader to App. [A] for a derivation, and simply report here the 
results. We remark that from now on we will work in the unitary gauge, where the Higgs doublet 



reads (0 , H/y2) . The effective coupling of the Higgs boson to two gluons reads (see App. \A 
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This coupling governs the rate of single Higgs production via gluon fusion, and its expression was 
already obtained in Refs. [10[|27|. The production rate normalized to the SM one is given by the 



square of the expression in square brackets in Eq. (11). On the other hand the effective coupling of 



two Higgs bosons to two gluons, which contributes to Higgs pair production via gluon fusion, has 
the following expression 
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(12) 



In terms of the effective Lagrangian in Eq. ((61) , the couplings read 
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(13) 



In the expression of the hhgg coupling in Eq. (13), the first term comes from the triangle top loop 



involving the tthh vertex, whereas the second is the contribution of top box diagrams, see Fig. |2j 
and kg and k 2g are parametrizing the contributions from integrated-out heavy particles. 



2 The bottom contribution is non-negligible, but can of course not be computed using the low-energy theorem, due 
to the smallness of the bottom quark mass. 
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3.2 Higgs interaction with photons 

Although the main focus of this paper is on gluon fusion, we give here the expression for the 
coupling of the Higgs boson to photons, as it is another loop process of crucial relevance for Higgs 
phenomenology at the LHC. This coupling receives contributions both from loops of heavy fermions 
and from the W boson loop. Application of the LET leads to the following effective Lagrangian [7| 



p 2 7 



-eff 



F^F^(J2Q 2 ftegm}(H) - -\ogm 2 w {H)) , (14) 



16tt 2 " J ' 4 

/ 



which is valid for nih < 2mw , 2m t , and where we have assumed that the heavy fermions transform 
in the fundamental representation of SU(3) C . Expanding around the VEV we obtain the /177 
interaction 



e 2 



-V- = ^2 F ^ h [Qt A i - I (-^logm 2 w (lI) ) J . (15) 



where we have assumed that all fermions have electric charge equal to that of the top quark[^J 
Qf = Qt = 2/3 , and A\ was defined in Eq. Q . By performing simple manipulations we obtain 
(see App. \M 
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(16) 



where we have replaced the LET approximation for the W loop with the full result encoded by the 
function 

J 1 {x) = F x (x), Fi(x) = 2 + 3x[l + (2- x)f(x)], f(x) = arcsin 2 (x- 1/2 ) , (17) 

which tends for large x to J 7 (oo) = 7 = 22/3 — 1/3 , where the first term comes from the transverse 
polarizations of the W and is precisely equal to the gauge contribution to the (3 function of the 
SU(2)l coupling, while the second term arises from the eaten Goldstone bosons. This replacement 



implies that the formal validity of Eq. (16) is extended to m;, < 2mj . The rescaling of the decay 
width T(/i — > 77)/r(/t — > 77)sm is obtained by comparing the square of the expression multiplying 
hF^ v F^ v in Eq. (16) in the two cases. In terms of the parameters of the effective Lagrangian in 
Eq. ^ the coupling reads 

2 7 

£hTY = ^2 F ^% ( 4( 3* c + kr r ~ aJ ^ m w/ m D) ■ ( 18 ) 



3.3 Single Higgs production via gluon fusion 

For a SM Higgs boson, the gluon fusion process 1 30 1 is the dominant production cross section at 
the LHC, see Refs. [3l] for reviews. At leading order (LO), the process proceeds via a top loop, 
with a subleading contribution from the bottom loop, see Fig. [TJ In composite Higgs models extra 

3 In all models considered in this paper, only top-like resonances contribute to the /177 coupling. The extension 
to heavy states with different electric charge is straightforward. 
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Figure 1: Generic diagram contributing to Higgs production in gluon fusion. In the SM case we have 
Q = t,b. In composite Higgs models with additional heavy fermionic resonances these add to the particles 
Q running in the loop. 



heavy, colored fermions with sizeable couplings to the Higgs boson are typically present, whose 
contributions to the gluon fusion process should in general be taken into account. It has been 
shown [9-11 32] that in explicit constructions the gg — )• h cross section (computed in the LET 



approximation) is insensitive to the details of the heavy fermion spectrum, i.e. it does not depend 
on the couplings and masses of composites, but only on the ratio v/ f, where / is the overall scale 
of the strong sector. This was found to be true both in models with partial compositeness and in 
Little Higgs theories. In fact, although the top Yukawa coupling receives a correction due to the 
mixing with resonances which depends on composite couplings, this contribution is exactly canceled 
by the loops of extra fermions, leading to a dependence of the gg —> h rate only on v/f. This also 
implies that the cross section can be obtained by simply multiplying the SM one by c 2 , where c is 
the rescaling of the top Yukawa coming only from the nonlinearity of the cr-model, and neglecting 
corrections due to fermionic resonances. 

Let us review how this cancellation arises. It is due to the fact that the determinant of the 
heavy fermion mass matrix takes the form 

det M 2 (H) = F(H/f) x P(A i ,M i ,/), (19) 

where F is a function satisfying F(0) = since the top becomes massless in the limit of unbroken 
electroweak symmetry, and P is a polynomial of the composite couplings Aj and masses Mj, but 



independent of H. It is then immediate to see that the hgg coupling in Eq. (11) does not depend 
on the masses and couplings of the resonances The origin of the factorization in Eq. (|19l) was 
explained in the context of partial compositeness in Ref. [11| , by means of a spurion analysis. There 
it was also pointed out that such a factorization can break down if the top mixes with more than one 
composite operator, leading to a dependence of the hgg vertex on composite couplings. Nevertheless 
in many explicit constructions, including Little Higgs models, the factorization in Eq. ( |19| ) takes 
place. Still, the independence of the hgg vertex on the composite couplings (collectively denoted 
by Aj) holds exactly only in the LET approximation, and corrections due to finite fermion mass 
effects are expected. We can estimate the residual dependence on the Aj due to finite fermion mass 
effects in a simple way. Assuming for simplicity the presence of only one top partner T, the mass 
eigenvalues can be written at 0(1/ f 2 ) as 



mt (H) = ^U-^ I pi , mT (H) = \ T f(l + a T 1 ^) , (20) 



The coefficient ch does not receive contributions from the heavy fermion sector. It is, however, generated by 



integrating out heavy scalars or vectors 27 
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where ax is a parameter dependent on the couplings Aj = {yt,^r} as a*r = 0(y% /A T )F] On the 
other hand, we can write = c y + 0(yf /X T ), where Cy is a constant arising from the pure 
nonlinearity of the a model. The low-energy theorem result for the hgg coupling reads, taking the 
effect of ch into account, 

= JL^^ f 1 - (4 ] - ^ + c -f\ o . (2i) 



48vr 2 ^ !)\ V y 2 J 

Notice that in the limit where T is heavy, corresponding to large \t, the effects of the heavy 
resonance on the hgg coupling vanish. In fact, ax goes to zero, whereas (£p — > Cy , implying that 
only the nonlinearity in the top Yukawa arising from the nonlinear a- model is relevant. 

By using the expression of the top Yukawa coupling (mt/v)(l — (cy + ch /%)£,) we can compute 
explicitly the top loop diagram, retaining the first subleading term in the 1/m 2 expansion. This is 



the leading correction to the LET coupling, given that my 3> nit- Thus Eq. (21) is improved to 



a z h 
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H 



(22) 



where we have used s = m?, and the ellipses stand for subleading corrections (including terms of 
order l/m T ). The independence of the LET hgg vertex of the composite couplings Aj is equivalent 
to the statement that Cy* — 2ay is a constant, — 2ax = c^j^] If this is the case then the 



dependence on the Aj of Eq. (22) is due to the last term, and we can estimate the sensitivity of the 



cross section to the Aj to be, for a light top partner At ~ yt, 

h £ ~ 0.06 C , (23) 



r^i 

a(gg -)■ h)sM 60 m 2 



where in the last equality we assumed rrih = 125 GeV. Thus corrections are expected to be very 
small even for large £. This estimate will be confirmed in Section [5j where the gg —> h cross 
section will be computed in MCHM5 retaining the full mass dependence. We note that in this 
model the cross section can be strongly suppressed compared to the SM, reaching ojosu = 1/3 
for a low compositeness scale £ = 0.25. However, the Higgs branching ratios into 77, WW and ZZ 
are enhanced compared to the SM, so MCHM5 can still be compatible with the excesses observed 
by ATLAS and CMS at rrih ~ 125 GeV even for values of £ as large as those considered in this 



paper 26 . 



3.4 Double Higgs production via gluon fusion 

Within the SM, double Higgs production via gluon fusion received interest mainly because it is 
sensitive to the trilinear Higgs self-coupling [33], see the first diagram in Fig. [2] In composite Higgs 
models, the process gg —> hh is affected essentially in two ways. First, the nonlinearity of the 
strong sector gives rise to a ffhh coupling (which vanishes in the SM) and thus to a genuinely 
new contribution to the amplitude, see the second diagram in Fig. [2j Second, one should take into 
account the effects of top partners, which include also new box diagrams involving off-diagonal 

5 If the quadratic divergence in the Higgs mass due to the top is cancelled by T, then the absence of an 0(H 2 ) 
term in Tr M 2 (H) = m 2 (H) + mr(H) implies ar = — J/t/(4Ay). See for example the explicit values of and clt in 



the Littlest Higgs model, reported in App. B Eq. (93 1 



e Notice that by using Eq. (|20| one finds det M 2 = y 2 \ 2 T f 2 H 2 (l - (c { y ] - 2a T )H 2 / f)/2. So if the factorization in 
Eq. (191 holds then Cj, — 2ar is a constant. 
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m rrii 



Figure 2: Generic diagrams contributing to double Higgs production via gluon fusion in composite Higgs 
models with rif novel fermionic resonances of mass (i = l,...,n/). The index j is introduced to indicate 
that different fermions can contribute to each box diagram. 



Yukawa coupling^] (shown in the second line of Fig. [2]). A first study of gg — > hh in composite 
Higgs models, neglecting top partners, was performed in Ref. [17| , where it was found that an 
enhancement of the cross section is possible due to the new tthh coupling (see also Ref. [16] for an 
earlier study in the context of Little Higgs models). For example, in MCHM5 with £ = 0.25, which 
corresponds to / ~ 500 GeV, the cross section was found to be about 3.6 times larger than in the 



SM. Recently, Ref. 18 performed a model-independent study of the process, making reference to 



the effective Lagrangian in Eq. ^ and thus again neglecting the effects of top partners, and found 
a large sensitivity of the cross section to the C2 coefficient parameterizing the tthh coupling. 

In this paper we include for the first time the effects of top partners in double Higgs production 
via gluon fusion. This is especially interesting in the light of the results of Refs. [4], where a light 
composite Higgs was shown to be tightly correlated with the presence of light top partners, as such 
light resonances can in principle affect the gg — > hh cross section in a sizeable way. Our analysis 
will confirm that this is indeed the case. 

We start by discussing the cross section in the LET approximation, which greatly simplifies the 
computation. In this limit, the amplitude is simply the sum of two diagrams, one with the effective 
hgg coupling followed by a trilinear Higgs coupling and the other involving the effective hhgg 
coupling. Adopting the SILH formalism, and recalling the expressions of the relevant Feynman 
rules 



2 \ d 



hgg : i ^8 ab {p^ - Pl ■ p 2 g^) 
hhgg : i g^(p&£ - Pi ■ P2 gT) [\ {{ai^W 
hhh : -*3=& [l + e(c6-|cH-|cr)] 



2 s 



jj logdetA^(F) 

/ H=v 



d log j 



jj ) log det M 2 (H) 



H=v 



(24) 



(where p± ; 2 denote the momenta of the incoming gluons), we can write the amplitude as 



-4let {gg ->■ hh) 



37TV 



5 ab (plp%- pi -P2 sDOunis), 



(25) 



7 Note that these Yukawa couplings only involve the top quark and its charge 2/3 composite partners. 
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where 



Clet(s) 



3m| 

9 

; — mi 



9 

; — mi 



1 ( d 

2 \dlogH 

Iff d 2 



logdet M 2 (H) 



d 



+ e (ce - 2c H 



i-£(c, 



c 6 + 2c# + — 



3c 



1 + ^{Cy + V " 3C 



(26) 



with s = (pi + P2) 2 denoting the partonic center-of-mass (c.m.) energy. It is immediate to check 
that the combinations c y — cq + 2ch + Or / 4 and c y + c r /4 are invariant under the reparameterization 
in Eqs. J2]) and ([3]). For completeness, we also give the result in terms of the effective Lagrangian 
in Eq. (l6|): 



Clet(*) 



The partonic cross section reads 



3m 



^2 (c + k g ) d 3 + 2c 2 - c 2 + h 



m 



29 ■ 



(27) 



!)!)- 



Vhh 



128(2vr) 3 9 



4m? 



C 



LET l 



(28) 



The hadronic cross section is obtained by convolution with the parton distribution function f g /p 
of the gluon in the proton, 



a 



r 1 r 1 d x 

/ dr —f g /p(x,Q) f g /p(T/x,Q)a gg ^ hh (Ts) 



(29) 



with the collider c.m. energy s related to s by s = ts. The renormalization scale \x and the 
factorization scale Q are chosen equal to the invariant mass of the Higgs boson pair, fi = Q = \f§. 
Throughout the paper, the parton distribution functions of MSTW2008 34 are employed. For 
£ — > 0, Eq. (26) correctly reproduces the SM result in the limit of large top mass 13, 14 , 



r SM m 



3m£ 

i — ml 



1 . 



(30) 



In the SM the mt — > 00 limit gives a cross section in agreement with the full result only within 20% 
for m h < 200 GeV (for m h = 125 GeV we find of^ = 14.6 fb and <rfj& = 17.9 fb) and moreover, it 
produces incorrect kinematic distributions, as noticed in Ref. [15]. Thus we expect the LET to be 
in general less accurate in gg — > hh than in single Higgs production. 



From Eq. (26) we read off that in models where the factorization Eq. (19) of det Ai 2 holds, the 
gg — > hh LET cross section is insensitive to composite couplings, due to a cancellation completely 
analogous to the one that we discussed for single Higgs production. In the left panel of Fig. [3] we 
show for mh = 125 GeV and a c.m. energy of 14 TeV the pp —> hh cross section normalized to 
the SM cross section (both were computed applying the LET) as a function of £ for some well- 
known models, in all of which the cancellation holds. We note that in MCHM5 the enhancement 
of the cross section is striking. This can be traced back to the behavior of the function Clet(s), 
which is proportional to the LET amplitude and is shown in the right panel of Fig. [3] for the three 
models under consideration. The enhancement for MCHM5 is evident. As pointed out for the 
first time in Ref. [17], where the gg — > hh process was studied in MCHM5 considering only Higgs 
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■ hh, 14 TeV, low-energy theorem 




800 1000 1200 1400 

m hh [GeV] 



Figure 3: (Left panel) The pp — > hh cross section for mh = 125 GeV at LHC14, computed using the LET, 
normalized to the SM cross section also computed in the m t — > oo limit. MHCM5 is discussed in detail 
in the text, whereas the gg — > hh amplitudes for MCHM4 and for the Littlest Higgs model are given in 
App.pl (Right panel) Square of the function C(m\ h ) 



, which was defined in Eq. ( 26 ) and is proportional to 



the LET gg - 
of m hh = Vs. 



hh amplitude squared, in the three models under consideration and in the SM, as a function 



nonlinearities (or equivalently in the limit of heavy top partners) but keeping the full dependence 
on mt, the dramatic increase of the gg —> hh cross section compared to the SM is mostly due to 
the presence of a new tihh coupling. The large enhancement of gg —> hh in MCHM5 is in contrast 
with the strong suppression in the same model of the single Higgs production cross section, which 
for £ = 0.25 equals 1/3 of the SM value (see Section [5]). 

By comparison with Ref. |17| we find that when fermionic resonances are above the cutoff, the 
LET underestimates the ratio (j M CHMh j a SM a |j OU t 30% . f or eX ample for £ = 0.25, application 
of the LET gives a cross section of 2.6 times the SM, whereas Ref. |17| found an enhancement factor 
of 3.6. This difference is due to the fact that in the former case nih <C mt is assumed, whereas in 
the latter the full mt dependence was retained. Notice that the best estimate of the cross section 
that can be obtained using the LET is a MCHM5 = (^let^ M V ct let) x a fuL because part of the 
corrections due to the finite top mass should cancel in the ratio of LET cross sections. In fact, 
in terms of cross sections the disagreement between the LET and the result obtained taking into 
account only Higgs nonlinearities is larger. For £ = 0.25 we obtain &^£ HM5 = 37 fb, whereas 
found (j MCHM5 = 64 fb, i.e. the difference is of order 50%. 
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Ref. 

In order to understand this behavior we investigated in more detail the validity of the LET both 
for single Higgs and double Higgs production. In single Higgs production the expansion parameter 
is m?/(4mf) and the series converges very quickly. In double Higgs production on the other hand, 
one needs to expand in s/(4m?) with s > 4m?, which is not small, so in general the expansion 
does not work as well as for the single Higgs case. In MCHM5, the validity of the expansion gets 
even worse. The reason why the LET is less accurate for MCHM5 than for the SM (where it 
underestimates the cross section by about 20%) is mainly the presence of the new triangle diagram 
containing the tihh coupling, which contrarily to the triangle diagram involving the virtual Higgs 
exchange does not vanish at large s. This is confirmed by taking into account the corrections of 
0(\jm\) to the LET result, which are reported in App. E.5 Compared to the SM we have an 



additional contribution ~ £ from the two-Higgs two-fermion coupling which goes like ~ s/mf, and 
which in contrast to the triangle diagram with virtual Higgs exchange is not suppressed by the 
Higgs propagator, see Eqs. (149) and (150). Therefore in MCHM5 for large £, where the coupling 
C2 is sizeable (see Table [I]), the corrections do not improve at all the LET result. For a model- 
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independent study of the gg — > hh process including only modifications to the top couplings, see 
Ref. [18]. 

While in the LET approximation the contributions of loops of top partners to the gg — > hh 
amplitude exactly cancel out with that coming from the modification of the top Yukawa due to 
mixing with resonances, the sensitivity to composite couplings of the full double Higgs production 
cross section (computed retaining all dependence on masses) is expected to be much larger than 
for gg — > h, where it was shown to be negligible. By direct computation in MCHM5, we will see 
in Section [6] that this is the case, i.e. the full gg — > hh cross section has a sizeable sensitivity to 
the details of the spectrum of the top partners. This effect is not captured by the simple LET 
result, which is completely determined by £. Therefore, while the low-energy theorem provides a 
useful tool to obtain a rough estimate of the cross section, a complete loop computation is needed 
to describe correctly the effects of top partners in the gg — > hh process. 



4 Composite Higgs model with extra fermionic resonances 



We consider a composite Higgs scenario with the symmetry group of the strongly-interacting sector 
given by 50(5), which is spontaneously broken down to 50(4) at the scale /. In order to cor- 
rectly reproduce the fermion charges an additional local symmetry U(l)x is introduced, leading 
to the symmetry breaking pattern 50(5) x C/(l)x/50(4) x U(l)x- This is the minimal realisa- 
tion including custodial symmetry and which implies four Goldstone bosons (GBs) transforming 
as a 4 of 50(4) ~ SU(2) L x SU{2) R . The SM electroweak group SU(2) L xU(1)y is embedded 
into 50(4) x U(l)x and the hypercharge Y is then given by Y = T| + X 
parameterized in terms of the field 



21 



22 



The GBs are 



S e 



n// 



n = -n/2T a /t a 



S = (0, 0, 0, 0, 1), (31) 

where T a (d = 1, ... ,4) are the generators of 50(5)/50(4) , and h a are the 4 real GBs. Using the 
explicit expressions of the 50(4) and 50(5)/50(4) generators = 1,...,5) 



rp a L,R 



rpa 

ij 



ubc/ x dec 



V2 



(5?6 



a^5 

3 



5* 5 5 



5fr) ± {8?5f - 
d = 1 4 



a L,R = 1,2,3 , a, b, c= 1,2, 3 



one obtains 



sm(h/f) 



h 



(h 1: h 2 , h 3 , h A , /icotan(/i//)) , h = h\ . 



(32) 
(33) 

(34) 



The two-derivative Lagrangian reads 
f 2 



-Ckii 



DfjTi = d M £ + ig WpXTt + ig' S M ST| . 



(35) 



By performing an 50(4) rotation, it is always possible to align the Higgs VEV to the /13 direction, 
thus identifying H = /13, where H is the Higgs field (with (H) 7^ 0). Then in the unitary gauge 



s 



co S (H(x)/f) 



S m(H(x)/f) 



1 



V 



EoC(z) 



sm(H(x)/f) cos(H(x)/f) ) 

;0, 0, sin(F//), 0, cos(H/f)), 



(36) 



(37) 
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and therefore 



C 



1 n 2 f 2 

kin= -d^H + ^f- sin 2 



2cos 2 9 w Z>lZ ^ 



(3i 



which fixes f 2 sin 2 ((tf ) //) = v 2 ~ (246 GeV) 2 . 

Fermionic resonances are described using the language of partial compositeness. We introduce 
vector-like fermions which have quantum numbers such that they can mix linearly with the SM 
fermions qi = (t^, &l) t and t R , and which at the same time have 'proto- Yukawa' interactions with 
the composite Higgs. We introduce composite fermions transforming as a complete 52/3 under 
SO (5) x U(l)x ■ This representation has the phenomenologically desirable feature that no tree- 
level corrections to the Z-b-b coupling arise, provided a discrete symmetry Plr exchanging the 
SU(2)l and SU(2) R factors is enforced (22l[35) . The composite multiplets can be written as 



1 



/ B-X h ' 3 \ 
-i{B + X 5 / 3 ) 

T + X 2 ' 3 
i(T - X 2 / 3 ) 
\ V2f j 



(39) 



Under SU(2) L x SU(2) R , a 5 of SO (5) decomposes as 5 ~ (2,2) (1,1) . The SU{2) L doublets 
Q = (T, B) T and X = (X 5 / 3 ,X 2 ^ 3 ) T form a bidoublet (2, 2) under SU(2) L x SU(2) R , while f is a 
singlet (1, 1) . The SM quantum numbers of the composite fields are summarized in Table ^1 Note 



field 


rp3 

1 L 




Y 


Q el = T 3 + Y 


T 


+1/2 


-1/2 


1/6 


+2/3 


B 


-1/2 


-1/2 


1/6 


-1/3 


X 5/3 


+1/2 


+1/2 


7/6 


+5/3 


X 2/3 


-1/2 


+1/2 


7/6 


+2/3 


f 








2/3 


+2/3 



Table 2: SM quantum numbers of the composite fermions in i\>. The last column denotes the electric charge. 

that Q has the same quantum numbers as the elementary doublet qL = (tL,^L) T , whereas T has 
the same quantum numbers as t R . The doublet X is peculiar of the 5 representation (it is absent 
in the most minimal case of the spinorial representation 4). Taking into account only one set of 
fermionic composites, the Lagrangian for the fermion sector then reads 



C f = iq L lfiq L +it R lfit R + ib R lpb R + i-ip L lft'ip L + i'ip R lfi'ip R 

- y/(^ L X T )(%) " MoMr + h.c. 

- A L q L Q R - A R f L t R + h.c. , 

where the covariant derivative acting on tp is given by 

= % - ig W^Ti - ig' B^(T R + X)] V , X = (2/3)l 5 . 
For later convenience, we also give the embedding of qL and t R in SO(5) vectors: 



(40) 



QL = V2 ^ L ~ ibL tL UL ' 



r/j = -J=(o V2t R ) T 



(41) 



(42) 
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Using these expressions the linear mixings can be rewritten as 



Cr. 



-A L Q L i) R - A r i/j l Tr + h.c. 



(43) 



From the Lagrangian in Eq. (40) we obtain the mass terms and the Yukawa couplings. The mass 
matrix reads 



Cr 



I 



\ 



ft 

ft 



\ 



J 



I 



\ 









A L 

yfs 2 

2 

yfsc 

V2 





yfs 2 

2 

yfsc 

V2 





yfsc 

V2 

yfsc 
M + yfc 2 



J 



( 



\ 



tR 

Tr 

Y 2/3 
X R 
Tr 



\ 



+ h.c. 



(44) 



where we have introduced the abbreviation s = sin((H) //) = v/f and analogously c = cos((H) / /). 
The diagonalization of the matrix, which mixes fundamental fields and composite states, is imme- 
diate before electroweak symmetry breaking, i.e. for v = 0. Then the mass terms become diagonal 
after the rotations 



Ql 

tR 

Tr 



cos <p£, sin <pL 

— sin 4>l cos tpi, 

cos 4>r sin 4 

- sin 6 R cos (p R 



(ft 
Ql 

t R 
Tr 



tan 4>l 



tan 4> R 



M 

Ar 
Mo + yf 



In this limit the top is massless, whereas the masses of the composite states are 



Mi 



Q 



Mo 

CL 



M X = M 



M* 



yf + M 

CR 



(45) 



(46) 



Electroweak symmetry breaking effects generate an additional mixing, which also involves and 
tR. Thus the top becomes massive due to its mixing with composite states. At the leading order 
in £ = v 2 / f 2 we have 



m t = y sin (p L sin (p R 



V2 



(47) 



Furthermore, the masses of the composite fermions in Eq. (46) get corrections of the order £. The 



Lagrangian Eq. (40), however, does not give rise to a mass for the bottom quark, because there is 



no composite in tp that has the right quantum numbers to mix with bR. Rather than introducing 
another fermionic multiplet (e.g., a 5_i/^) to solve this issue, we introduce a small breaking of the 
partial compositeness pattern, namely a Yukawa coupling of the Higgs to elementary states 



Cb 



-X b q L Hb R + h.c. , 



(48) 



where H is the Higgs doublet. We will, however, neglect the small effects proportional to . 

Expanding the proto- Yukawa term up to second order in the physical Higgs h we obtain the 
Higgs Yukawa couplings and the two-Higgs two-fermion couplings. The Yukawa coupling part of 
the Lagrangian reads 



C hff =yh 



( \ 

T L 



fo 




V° 





sc 
sc 

l-2s 2 
V2 





SC 

sc 

l~2s 2 
V2 



o \ 

l-2s 2 
l-2s 2 

-2sc J 



( ** \ 

Tr 

Y 2/3 
X R 

\f R J 



+ h.c. 



(49) 



G 



hff 
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For the two-Higgs two-fermion interaction Lagrangian we find 



£hhff 



-A" 
V 



( tL \ 






Vo 





1 - 2s 2 
1 - 2 s 2 
-2y/2sc 




-2 s 2 -2\/2sc 
-2V2sc 
-2(1 -2s 2 ) J 



1 

1 - 2s 2 
-2y/2sc 



( \ 

Tr 
Y 2/3 
X R 
\f R J 



+ h.c. 



(50) 



G 



hhff 



After the rotation into the mass eigenstate basis, the two matrices G hf f and G hhf j yield the single 
and double Higgs couplings to fermions, respectively, which will be needed for the calculation of 
the single and double composite Higgs production cross sections through gluon fusion. 

4.1 Constraints from electroweak precision data and flavor physics 

The strongest experimental constraints on composite Higgs models still come from the electroweak 
precision measurements at the Z pole mass at LEP. A convenient description of LEP precision data 
is given in terms of the parameters e±, 62, £3 and [36] - |38| . These parameters are on the one hand 
measured with high precision [39j , and on the other hand can easily be computed theoretically. In 
addition to the SM contribution present in the decoupling limit / — > 00, the MCHM5 contributes 
to the e parameters through three different effects. The first beyond the SM (BSM) effect arises 
from the modified coupling of the Higgs to W and Z gauge bosons, which induces a logarithmically 
divergent contribution to the oblique parameters S and T, or equivalently to ei and £3. The 
contribution is cut-off by the mass m p of the first composite vector resonance |40| , 



Ae 1 * 




16-7T cos 2 6w 



Ae 1 * 



a(M z 



487r sin 2 



-Clog 




(51) 



The second effect is the direct contribution of the vector p and axial-vector a resonances to the S 
parameter, which in the MCHM is found to be (see Ref. |41j) 



in 



w 



mt 



1 + 



1.36 



m 



mt 



mt 



(52) 



In the second equality, we have used the relation m a /m p = 5/3, obtained in the five-dimensional 
realisation of the model 21 . The third and last contribution to electroweak precision parameters 



comes from the top partners at one loop, giving contributions both to the T parameter and the Z-b-b 
vertex, i.e. respectively to e\ and (24,42-44 . Computing the precise value of these contributions 



requires the numerical diagonalisation of the mass matrix of the top quark and its partners, which 
depends on the parameters A^, An, Mq, y and /. The requirement that the top mass matches the 
measured value mt = 173.3 GeV allows, however, to express the corrections to ei and e& in terms 
of four dimensionless parameters, 



Ae[ 



fermions 



fl (£,(/>L,<t>R,R) , 



a fermions 

Ae J b 



fb [£,<t>L,<t>R,R) , 



(53) 



where £, 4>l, 4>r are defined above and R = (Mo + yf) /Mq. The function fi is computed exactly 
at one loop, while for fb only the longitudinal polarisations of the gauge bosons are taken into 
account in the loop. The values obtained in this way are consistent with the full one-loop result 
of Ref. 44 . The agreement of the model with experimental data is then assessed through a x 2 



test, described in detail in App. O The latest electroweak precision data are used, including the 
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2012 update of the W mass. Fixing the Higgs mass to rrih = 125 GeV, the model is completely 
determined by the five parameters £, (j>z, <j)R, R and m p . Over the latter four a scan is performed for 
£ fixed to two representative values, namely 0.1 and 0.25. The results are displayed in Fig.|4]for the 
left-handed compositeness angle 4>l versus the mass of the lightest top partner. Note that the value 
of R is bounded by the requirement that y < An and that we impose furthermore the constraint 
\Vtb\ > 0.77 |45|. In Fig. [5j the whole spectrum of composite fermions is shown for a sample of 
parameters passing the EWPT. The green points correspond to B, the red ones to X 5 / 3 , and the 
blue points for each set of parameters denote to the top partners T, X 2 ! 3 and T, which cannot be 
properly distinguished one from another once the rotation in the physical basis is performed. At 
leading order in v/f, composite fermions within an electroweak doublet have, however, the same 
mass, so that the green points describe effectively the mass of the (T, B) doublet and the red ones 
the mass of the (X 5//3 , X 2 / 3 ) doublet. The blue points far from the red and green regions in Fig. [5] 
can therefore be interpreted as singlets T. 

There are two regions of the parameter space compatible with EWPT where in addition at least 
one of the top partners is light, as generically needed in order to have a light enough Higgs. The 
first region corresponds to low values of the top compositeness angle 4>l, where the lightest top 
partner is typically the singlet T. In this case the fermion bidoublet is always heavier than 1.5 TeV 
(see Fig. [5| and decouples. Note that the right-handed top must then be very composite in order to 
yield the correct top Yukawa coupling y t . The second region corresponds to large values of sinc^x,, 
for which the top-bottom doublet becomes fully composite. In this second region, the 'custodian' 
doublet X is very light, having a mass well below a TeV. Since the X doublet contains an exotic 
charge 5/3 fermion (which turns out to be the lightest new fermion for large sin^^), this region is 
very sensitive to direct collider constraints, as will be discussed in the next section. An intermediate 
region with moderate values of sin</>£ is also allowed by precision data, although all new fermions 
are rather heavy, above 1 TeV, so it could be difficult to obtain a light enough Higgs in this region. 
A thorough discussion of the implications of top compositeness is contained in Ref. [46]. A comment 
is however in order: in Ref. |46|, when mixes with a (1,1) of SU{2)l x SU(2)r (which is the 
case in MCHM5) a highly composite right-handed top is not accompanied by anomalously light 
fermionic resonances, because there is no custodial partner of the state mixing with tR, the latter 
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Figure 4: A sample of parameters passing the x 2_ test of electroweak precision observables, displaying the 
compositeness of the left-handed top versus the mass of the lightest top partner, for £ = 0.25 (left) and £ = 0.1 
(right). The points in light gray do not pass the direct collider constraints. Points in orange/medium gray 
pass the present const raint s but will be tested by the LHC running at 8 TeV with an integrated luminosity 



of 15 fb , see Section 4.2 
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Figure 5: Physical mass spectrum of the composite fermions for a sample of points passing the electroweak 
precision tests, as a function of the left-handed top compositeness for £ = 0.25 (left) and £ = 0.1 (right). The 
blue/dark gray points are top-like fermions (charge +2/3), the green/fair gray points bottom-like (charge 
— 1/3), and the red/medium gray ones correspond to the exotic X (charge +5/3). Light gray points are 
excluded by present collider constraints, see Section [42} 



being T in MCHM5. Indeed one can write from Eq. (46) Mf = Ar/ sin^/j — > Ar when sin c/)r — > 1, 
implying that T is not necessarily light in presence of a strongly composite tR. However, EWPT 
select a light T for large sin^, as discussed for example in Ref. [431. 

Note finally that the constraints on the parameter space from electroweak precision data can 
be significantly relaxed by extending the fermion sector of the model [44] . 

Additional constraints on the model come from flavor physics. Composite Higgs models generi- 
cally give rise to four- fermion operators which contribute to flavor-changing processes and to electric 
dipole moments. Low values of the compositeness scale / as considered in this paper are only al- 



lowed if the strong sector is flavor-symmetric, hence effectively implementing MFV 47 . In this case 
both flavor-changing processes and electric dipole moments are inhibited, but the MFV assumption 
requires a large degree of compositeness also for light quarks, which are therefore sizably coupled 
to the strong sector resonances, leading to a different phenomenology. 

Experimental constraints can be described in an effective formalism, in which four-fermion 
operators arise after integrating out the vector resonances. The most relevant operator is 

' ,Z (sin <^) 4 (q L ^t a q L ) (q L7 ,t a q L ) , (54) 



4m^ 



where t a are the generators of SU(3) C , which imposes a constraint on the size of the mixing angle 
4>l- From the most recent experimental dijet angular distribution [48], the bound is 

\(sin^ L ) 4 < — ^— - 2 ( Sin(/)i )2<_f_ (55) 



4m2 v YJ ^> ~ (7.5 TeV) 2 ~ 1.5 TeV ' 

or equivalently sin</>£ < 0.6 for £ = 0.25 and sin^ < 0.7 for £ = 0.1. Similar bounds apply to 
the compositeness of right-handed quarks. An even more stringent constraint comes then from the 
compositeness of the left-handed bottom quark, which affects the mass splitting between neutral 
B mesons, yielding the bound [8] 

£ (sin^) 4 < 2-HT 3 . (56) 

Flavor constraints are therefore strongly disfavoring the region of the parameter space where 
sin <f)L ~ 1, corresponding to a light doublet of top partners. They rely nevertheless on a number 
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of estimates, and we will therefore not consider this region as fully excluded. 



4.2 Constraints from searches for heavy fermions 

Expanding the composite Higgs coupling to fermions, we obtain the leading interactions between 
one heavy fermion and two SM particles, which mediate the decay of the heavy states, 



£ v t/*L<'-u ( b L ir + -^={h ~ r ''T 



°)J Tr + V^2/3L [h + i vr ) t R + h.c. 



+ -^s R c L T L (h - i7T°) t R + ys R c L B L t R ir - ys R X 5/3L t R ir + + h.c. , (57) 

where we have already performed the rotations oc 4>l,r ■ From the Goldstone equivalence theorem 
then follow the leading order branching ratios (in the limit S> mz, mh) 

BR(T -»• Wb) = i , BR(T -»• Zt) = BR(T ->• hi) = ^ ; 
BR(X 2/3 -»• Zt) = BR(X 2/3 -»• Tit) = - , BR(X 5/3 -»• Wt) = 1 ; 
BR(T -> Zt) = BR(T -»■ /tt) = - , BR(£ Wi) = 1 . 
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Figure 6: Branching ratios of the lightest top partner into W + b (upper left), Zt (upper right) and ht (lower) 
as a function of its mass for £ = 0.25. Points for which the singlet (doublet) is the lightest top partner are 
shown in pink/medium gray (purple/dark gray). 
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However, in our analysis of the electroweak and collider constraints we will keep all orders in £, by 
performing a full numerical diagonalization of the mass matrix in the top sector and computing 
the couplings of the mass eigenstates to gauge bosons and to the Higgs boson in the unitary gauge. 
Complete formulae for the partial decay widths of heavy fermions into SM fields are given in 
App. |Dj As we will see, for the relatively large values of £ that we consider, sizeable corrections to 
the leading approximations listed above arise. In Fig. [6]we show the branching ratios of the lightest 




Figure 7: Cross sections for QCD pair production of heavy fermions at approximate NNLO, at the Tevatron 



(dashed), at LHC7 (solid) and at LHC8 (dot-dashed). The cross sections were computed using HATHOR 49 
and MSTW2008 PDFs. 

top partner as a function of its mass for £ = 0.25. The pink points are the branching ratios in case 
the lightest top partner passing the electroweak prescision tests is the singlet. The purple points 
correspond to the lightest partner in accordance with the EWPT being the doublet. Compared to 
the approximate formulae the branching ratios into ht are a bit enhanced while the ones into Zt 
are somewhat reduced. 

In certain regions of the parameter space, some fermionic resonances can be very light, thus 
rendering constraints from direct searches for heavy fermions at the LHC and Tevatron relevant. 
The experimental collaborations have performed several searches for pair-produced heavy fermions, 
with subsequent decay into the final states WbWb, ZtZt, WtWt . Since pair-production of the 
heavy fermions is a QCD process, the cross section cr(pp, pp — > tpif)) , with ifi being a generic heavy 
fermion, only depends on M^u . The constraint from e.g. a search for ipip — > WbWb at the LHC 
will read 

<?qcd{pp^W^) xBR(V> ^Wbf <a exp , (58) 

where a exp is the upper bound on the cross section, as given by the experiment for each value 
of the resonance mass. The QCD pair-production cross sections were obtained at approximate 
next-to- next-to-leading order (NNLO) [49], and are shown in Fig. [7j We remark that also single 
production of heavy fermions can give complementary, relevant constraints (see e.g. the fourth 
among Refs. |5j for a detailed analysis), however no such search has been published by ATLAS and 
CMS yet. Therefore, we do not discuss single production. 

Note that the branching ratios are non-trivial only for the top partners T, X 2 / 3 and T, whereas 
B and X 5 / 3 decay with unity branching ratio into tW T . We summarize in Table [3] all the searches 
for pair-produced heavy fermions that we included in our analysis. The analyses of tWtW final 
states, although intended by the experiments to be searches for heavy charge —1/3 quarks such as 
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cxp. 


search 


L [fb- 1 ] 


range in My, [GeV] 


ref. 


CMS 


50 




WbWb (1 lepton) 


4.7 


[400, 625] 


CMS-PAS-EXO- 11-099 








WbWb (2 leptons) 


5.0 


[350, 600] 


arXiv:1203.5410 






wtwt 


1.14 


[350, 550] 


CMS-PAS-EXO- 1 1-036 






WtWt 


4.9 


[450, 650] 


arXiv: 1204. 1088 






ztzt 


1.14 


[250, 550] 


arXiv:1109.4985 


ATLAS 


51 


WbWb 


1.04 


[250, 500] 


arXiv: 1202.3076 






WqWq 


1.04 


[300, 500] 


arXiv: 1202.3389 






WtWt (1 lepton) 


1.04 


[300, 600] 


arXiv: 1202.6540 






WtWt (2 leptons) 


1.04 


[300, 600] 


arXiv:1202.5520 


CDF 


52 




WbWb 


5.6 


[180,500] 


arXiv: 1107.3875 








WtWt 


4.8 


[260,425] 


arXiv:1101.5728 



Table 3: List of experimental searches for pair-produced heavy fermions that we included in our analysis of 
collider constraints. 

the B, apply straightforwardly also to the X 5 ^ 3 , which decays into the same final state^} 

The region of the parameter space corresponding to sin (f>L ~ 1 is the most constrained by direct 
searches. The lightness of the X 5 ^ 3 fermion in this case (see Fig. [5J is prohibited by both Tevatron 
and LHC searches in WtWt final states. For a lower degree of compositeness of the left-handed 
top quark, the lightest top partner is the singlet T, which decays in all three final states Wb, Zt 
and ht. The Tevatron only has enough sensitivity to exclude top partners below 300 GeV, while 
the most stringent LHC constraints {i.e. those based on the full 2011 luminosity) start only at 
350 GeV. This leaves a region of the parameter space in the range rrif G [300, 350] GeV which is 
not directly excluded by present constraints, see Fig. |4} 

In addition to the present exclusion limits, we show in Fig. [4] an estimate of the reach of the LHC 
in 2012. The increase in energy enhances significantly the production cross section of heavy fermion 
pairs (see FigjT]). On the other hand, the present exclusion limits quoted by ATLAS and CMS will 
be modified due to the changes in the background and to the additional integrated luminosity. 
Backgrounds in searches for top partners are dominated by top quark pair production, which is 
increased by 42% when going from 7 to 8 TeV cm. energy at the LHC. The search strategy relies 
on a cut on the tt invariant mass, whose distribution is not significantly affected by the increase in 
energy, as explicitly checked using MadGraph 5 [53] . The upper limit on the top partner production 
cross section is therefore softened in the Gaussian approximation by a factor vL42 = 1.19. The 
total luminosity of 15 fb -1 expected to be attained in 2012 is nevertheless tightening the limit on 
the cross section, lowering it by a square root factor of the luminosity in every channel. More refined 
searches after the LHC upgrade to 14 TeV will be needed in order to explore the full parameter 
space (5]. 

5 Single Higgs production in MCHM5 

The cross section for single Higgs production in MCHM5 can be readily derived by noting that 
we can directly apply Eq. ^ since the Higgs kinetic term is canonically normalized. We therefore 

8 Note that the decay products of BB and X s ^ 3 X 5 / 3 would have different spatial configurations. For example, 
same-sign leptons necessarily stem either from X 5 ^ 3 or from its antiparticle, while in the case of BB production each 
of the same-sign leptons arises from a different heavy particle. However, since in the current searches only basic cuts 
on single objects are applied, this kind of kinematic differences is expected to give negligible effects on the exclusion 
limits. 
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only need the determinant of the mass matrix of top-like fermions in Eq. ( 44 ) , which reads 



8 cos z 6 1 



b L sin 2 (j> R 



8 COS 2 (pL COS 2 4>R 



(M + y /) 2 sin 2 



(59) 



and has the form of Eq. (19) 
sm 2 ((H)/f) = 0, and 



Thus we obtain A\ = (2/v)(l — 2£)/- v /l — f (where we have used 



a(pp h) 



a(pp -> h)sM W 1 - I 



1 - 2£ 



(60) 



which is valid to all orders in £. Equation (60) is independent of the details of the fermion spectrum. 



While this holds exactly only in the low-energy theorem limit, as discussed in Section 3.3 we expect 



that retaining the full mass dependence will give corrections to the cross section at most of a few 
percent. This is confirmed by a full computation of the cross section in which the dependence on 
the fermion masses is retained, as shown in Fig. [8| The figure shows the cross section of single 
Higgs production through gluon fusion including new fermionic resonances, normalized to the SM 
cross section computed with finite mj, as a function of the mass of the lightest resonance^] Note 
that the QCD if-factort 10 cancel out under the assumption that the higher order corrections are 

(see Ref. 32 ). A parameter scan has been performed, selecting only 



n 



the same in both case; 

points that satisfy EWPT. Among these, points that satisfy all current collider bounds are shown 
in green, whereas gray points are already excluded. In orange we show points currently allowed, 
but that will be excluded by searches for heavy fermions at the end of the LHC8 run if no excess 



is observed. The agreement with the prediction of the low-energy theorem in Eq. (60), shown as 



a black line, confirms that the cross section is to an excellent approximation independent of the 
details of the spectrum, and is fixed only by £. The sensitivity to the composite couplings is at 
most 2% x osm for light top partners, in agreement with our previous estimate, and vanishes for 
heavy partners. We conclude that for single Higgs production the LET provides a very accurate 
cross section for any spectrum of the extra fermions. 



Finally we remark that the result in Eq. ( 60 ) coincides with the one obtained considering only 
the Higgs nonlinearities, i.e. rescaling the SM cross section by c 2 , where c is the correction to the 
top Yukawa in the limit where fermionic resonances are heavy and thus their effects negligible (see 
Table [I}. This is a consequence of the cancellation discussed above. 



5.1 Effect of non-minimal operators 

We can add to the minimal partial compositeness Lagrangian in Eq. ( |40[ ) the following operators 

A£ = i y' L $ L Y?)U){^ L ) + i ^(V^£ T )$(&/^) , (61) 

where the covariant derivative reads = — ig'X (X = 2/3). The most convenient way to 
discuss these operators is to perform the following field redefinition, 



L,R 



C(x) T 1pL,R ■ 



(62) 
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9 We have compared our results in the SM limit to the ones obtained with HIGLU 
10 The i^-factor is defined as the ratio of the higher-order cross section to the leading 



j order cross section. 

11 This assumption is valid only at next-to-leading order (NLO) QCD. At NNLO QCD different mass scales play a 
role. Furthermore, the correct matching of the strong coupling constant in the effective and the full theory has to be 
performed. In Ref. 32 , however, it was shown that for parameters similar to ours the differences in the /("-factors 



are of the order of a few percent only, so that the SM NNLO /("-factor can safely be applied also to the MCHM5 case. 
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Figure 8: The MCHM5 cross section for single Higgs boson production through gluon fusion (including the 
exact dependence on top and heavy fermion masses) , normalised to the SM cross section (computed retaining 
the mt dependence), as a function of the mass of the lightest fermion resonance miightest for TUh = 125 GeV. 
The compositeness parameter has been chosen £ = 0.25. Green/dark gray points are allowed, gray points 
are excluded by current collider constraints, whereas orange/fair gray points will be tested by LHC8 in 2012. 
For comparison, the cross section ratio computed with the LET, Eq. (l60|), is shown as a black line. 



where C( x ) was defined in Eq. (36). Upon this transformation, the Lagrangian reads (omitting 
kinetic terms of elementary fields, and gauge interactions) 



- y/(^£rJ)(£o^R) " M^ L ^ R + h.c. 
+ iy' L (^ L ^)lpCE ^ L ) + (L -> R) 

T I 



AlQlC [x)it> R - A R i; L ((x) Tr + h.c. 



(63) 



Thus we need to rescale the singlet T to make it canonically normalized, Tl,r — > Ti^/wl + y' L R . 
Let us now focus on how the amplitude for gg — > h is modified by the new operators. It is easy 



to verify that the Higgs derivative interactions contained in Eq. (63) do not contribute to the 
amplitude for single Higgs production, because they are antisymmetric in the fermion fields 
Therefore we can simply apply the low-energy theorem. From the fermion mass matrix, which 
i*6cicl.s 

A T l+cos(H//) a cos(H//)-l a sin(ff//) 1 \ 



/ 



M 







>/2, 



V2 
cos(ff//) 





Mo 







Ar 






^" 2 



M 




V2 y/T+jjr 



Mp+yf 

v /T +^v /T +^ / 



(64) 



we obtain 



1(1 + j/ L )(l + j/ R ) 



det M 2 (H) 



sm 



2H\ 
TJ 



(65) 



which implies that the amplitude for gg — > h is not sensitive to the value of y' L R , see Eq. (19). 



On the other hand, the Higgs derivative interactions in Eq. (63) contribute in general to the pair 
production process, because they enter box diagrams. Therefore the cross section for gg — > hh 
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will be sensitive to y' L d . In the following section, however, we consider the minimal Lagrangian, 
setting y' L R = . 

Finally we comment about the contribution of the exotic state X 5 ^ 3 to the amplitudes for 



gg — > h, hh. In the field basis of Eq. (63) the Higgs appears only in elementary/composite mixing 
terms and in derivative interactions (thus showing manifestly its pseudo-GB nature). Since 

X 5/3 

does not mix with any elementary field, there is no contribution to the amplitudes from the mixing 
terms. On the other hand, it is easy to check explicitly that Higgs derivative interactions do not 
involve X 5 / 3 . We conclude that the exotic state does not contribute at all to the amplitudes for 
single and double Higgs production in gluon fusion. 

6 Double Higgs production in MCHM5 

In this section we discuss the cross section for pp — > hh first in the LET approximation, and 
subsequently retaining the full dependence on the masses of the fermions running in the loops. 

6.1 LET cross section 



From the determinant of the fermion mass matrix in Eq. ( 59 ) we can compute A 2 = (-2/V)/(l-£), 
which determines the hhgg coupling via fermion loops. This, together with the form of A\ previously 
derived and with the expression of the /i 3 coupling given in Table [TJ allows us to write down the 
amplitude for gg — > hh at all orders in £: 



LET 3m 2 h ( 1 - 2£ \ 2 



Thus analogously to single Higgs production, the LET cross section for Higgs pair production is 
insensitive to the details of the heavy fermion spectrum, and is fixed only by £. The corresponding 
pp —> hh cross section at LHC14, normalized to the SM cross section (also computed in the infinite 
rat limit) was shown as a function of £ in Fig. [3] (left). 

6.2 Enhancement of the cross section 

We have seen that for small values of £ single Higgs production in the MCHM5 is suppressed 
compared to the SM while double Higgs production is enhanced. The behavior of single Higgs 



production becomes clear from the LET result given in Eq. (60). In double Higgs production the 
Higgs pair is either produced through Higgs bosons coupling to the gluons through triangle loops 
or through boxes. In the former case, in the SM we only have a diagram with a Higgs subsequently 
decaying into two Higgs bosons, while in composite Higgs models there is an additional triangle 
diagram due to the two-Higgs two-fermion coupling. In the amplitude for Higgs pair production 
the parts coming from the triangle containing the triple Higgs coupling and from the box diagrams 
interfere destructively. In the MCHM5 amplitude where these two contributions are modified by 
((1 - 2£) / 'y/Y^l) 2 , the additional dia gram with the two-Higgs two-fermion coupling proportional 
to £ can hence have order one effects so that it governs the total cross section. This can be inferred 
from Fig. [9] which shows the double Higgs production MCHM5 cross section normalised to the SM 
as a function of £ for three different approximations. The red line has been obtained in the limit 
of heavy top partners keeping the full top quark mass dependence, the blue line is the LET result, 
and the black line, finally, is obtained by explicitly setting the two-Higgs two-fermion coupling to 
zero. In this case the cross section ratio is given by ((1 — 2£)/ y/1 — £) both for the LET and for 
the approximation where the top quark mass dependence has been kept. The dotted lines in the 



25 




Figure 9: The cross section for double Higgs production in MCHM5 normalised to the SM as a function of 
£ for three different approximations. Red: In the limit of heavy top partners keeping the full top quark mass 
dependence. Blue: LET. Black: setting the two-Higgs two-fermion coupling to zero. The red/blue dotted 
lines show the same as the red/blue full lines after application of an invariant mass cut of m^t > 600 GcV. 



figure have been obtained by applying an invariant mass cut of m^h > 600 GeV. After application 
of the cut the discrepancy in the cross section results for the two approximations becomes even 
worse, see also the discussion in Section 16.4 



6.3 Full 1-loop cross section 

In the triangle diagrams which contribute to double Higgs production the gluons couple to the total 
spin S z = along the z-axis, whereas the box diagrams involve S z = and S z = 2 couplings. The 
amplitude for the process can hence be expressed in terms of independent form factors F/\, Fq, 
-Fn,5 associated with spin and Gn,5 associated with spin 2. The total partonic cross section 
is given by 



1 



Vgg^hh 



1024(2vr) 3 s 2 



+ 



dt 



4 4 



i=l j=l 
4 



(67) 



with the integration limits 



(68) 



where s denotes the partonic cm. energy. The triangle and box form factors as well as their 
expansion in 1/m 2 are given in App. |eJ The various couplings appearing in Eq. (67) are also 
defined there. We have explicitely verified that in the SM limit our result agrees with Ref. [14] . 
The hadronic cross section is obtained by convolution with the parton distribution function of the 



gluon in the proton, see Eq. (29). 
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Figure 10: The cross section for double Higgs production through gluon fusion normalised to the SM as 
function of the mass of the lightest resonance of the heavy top sector, for to/, = 125 GeV. The compositeness 
parameter has been chosen £ = 0.25. Green/dark gray (gray) dots denote points which pass (do not pass) 
all current constraints, whereas orange/fair gray dots correspond to points that will be tested by LHC8. 
The left panel shows points for which X 2 / 3 is the lightest top partner (as a consequence of being largely 
composite), whereas for points in the right panel the lightest top partner is typically the singlet T. The black 
solid (dashed) line corresponds to the result in the limit of heavy top partners keeping the full top mass 
dependence (to the LET result as in Fig. [5]) . The expected number of events in the hh —> bbjj final state 
after all cuts at LHC14 with L = 300 fb -1 is also shown, along with the 3cr evidence threshold (dot-dashed 
line), see text for details. 



6.4 Numerical analysis 

For the numerical analyis we have performed a scan in the parameter set (<f>L, <Pr, R) and retained 
only the points which fulfill the constraints from EWPT. By this we mean that there exists a value 
of m p G [1.5 TeV, 4m f] such that the configuration (/, 0£, <f>R, R, m p ) passes EWPT at 99% CL. For 



this set of points we show in Fig. 10 for nih = 125 GeV and £ = 0.25 the double Higgs production 
cross section normalised to the SM as a function of the lightest top partner mass. The dependence 
on the masses of the loop particles has been fully taken into account. The black solid line shows the 
result in the limit of heavy partners, keeping only the top contribution (with full mass dependence) 
in the loop, while the black dashed line corresponds to the LET result in Fig. [3j The green (gray) 
dots are points which pass (do not pass) the current constraints from Tevatron and LHC data, 
whereas orange points will be tested by LHC8. 

Some comments are in order. First of all, we find a sizeable dependence of the cross section 
on the spectrum of the heavy fermions with 2.7 < ct/ctsm ^$ 3.7. We recall that both the LET 
cross section and the cross section in the limit of heavy partners only depend on £, The LET 
approximation, however, severely underestimates the ratio cr/asM, and this effect is even worse if 
we refer directly to the cross section, since we are consistently normalizing the LET cross section 
for MCHM5 to crsM( m t — > 00), which is ~ 20% smaller than the full result. On the other hand, the 
result in the limit of heavy partners, while keeping the full top mass dependence [l7| , overestimates 
the cross section in the region mu g htest i$ 1 TeV, which is compatible with a Higgs as light as 
125 GeV. For large values of the partner masses, of course, the cross section tends to the value 
obtained including only top loops (with top couplings following the 'trigonometric' rescalings given 
in Table [Tf. 

It should be noted that we have not taken into account higher-order QCD corrections. They 
have been calculated at NLO for SM and MSSM Higgs pair production in Ref. |55| in the heavy 
top mass limit. However, they cannot be taken over here as we have the additional diagram with 
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the two-Higgs two-fermion coupling and more seriously box diagrams with different loop particle 
masses. For heavy loop particle masses we do not expect the corrections to be too different from 
the SM case, so that they approximately cancel out in the ratio of the two cross sections. 
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Figure 11: Cross section for double Higgs production through gluon fusion normalised to the SM as function 
of the mass of the lightest resonance of the heavy top sector, for £ = 0.1 and = 125 GeV. Points are split 
in the two panels depending on the degree of compositeness of tz- The meaning of the dots and lines is the 
same as in Fig. 10 except that we assumed an integrated luminosity L = 3ab _1 at LHC14. The dot-dashed 



line corresponds to the 5(7 discovery threshold, see text for details. 



In Fig. 11 we show the corresponding results for a lower value of £ = 0.1, which corresponds 
to / ~ 800 GeV. Due to the larger value of /, the cross section is less enhanced compared to the 
SM. Similarly to the case £ = 0.25 the LET underestimates the cross section, although in a less 
dramatic way than in the previous case. 

To estimate the reach of the 14 TeV run of the LHC on Higgs pair production, we focus on the 
hnal state hh —> 6^77, which was shown to be the most promising for a light Higgs boson 15 56-58 . 
Reference [58] found that assuming a luminosity L = 600fb _1 , 6 signal events could be obtained 
after all cuts, with a background of 11 events. We estimate the expected number of signal events 
for MCHM5 by computing a(pp — > hh) x BR(hh — > 6677) for each point in the parameter space 
(taking into account the QCD production iT-factoi 12 of 1.9 and the non-standard Higgs branching 



ratios) and multiplying it times the acceptance for all cuts as computed in Ref. |58| for the SM. 
This rough approximation cannot of course replace a full analysis of the effects of cuts in the 
MCHM5 case. The latter, however, goes beyond the scope of this work. We therefore apply the 
simplified procedure for an illustratory purpose. We also quote the number of events needed for a 
3 (5)cr evidence with L = 300 (3000) fb -1 , based on the background estimate of Ref. [58] with the 
requirement of one 6-tagged jet. Notice that this is likely to be conservative, because the analysis 
of reducible backgrounds (whose sum is larger than the irreducible 6677) performed in Ref. 58 



made use of efficiencies and misidentification probabilities, in particular for 6-jets, that have since 
then been improved by ATLAS and CMS. We find that a 3a excess can be obtained already with 
300 fb -1 if £ = 0.25, except perhaps in some regions of the parameter space with a very light top 
partner. A 5a discovery would be possible at the LHC luminosity upgrade for a more moderate 
value £ = 0.1. 

We note that in the recent Ref. [18] two 6-tagged jets were required, and the efficiency and 
rejection probabilities for 6-tagging updated to current values. However, since we are only inter- 

12 As stated above the SM QCD corrections to double Higgs production cannot be translated trivially to the 
composite Higgs case. Assuming the top partners to be heavy we expect, however, the error not to be too large by 
applying the SM K-factor to MCHM5 double Higgs production. Concerning the diagram involving the two-Higgs 
two-fermion coupling we explicitly verified that it hardly changes the QCD corrections compared to the SM ones. 
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Figure 12: The cross section for double Higgs production through gluon fusion normalised to the SM for 
to/j = 125 GeV, as function of the mass of the lightest resonance of the heavy top sector. The compositeness 
parameter has been chosen £ = 0.25. Green/dark gray (gray) dots denote points which pass (do not pass) all 
current constraints, whereas orange/fair gray dots correspond to points that will be tested by LHC8. The 
left panel shows points for which X 2 / 3 is the lightest top partner, whereas for the points in the right panel 
the lightest top partner is typically the singlet T. The black solid line corresponds to the result obtained 
considering only Higgs nonlinearities, i.e. in the limit of heavy top partners and keeping the full top mass 
dependence. 



ested in a rough estimate of the LHC reach, we conservatively adopt the numbers of Ref. |58| . 
Furthermore, a realistic analysis of the instrumental backgrounds relevant to 6^77 would require a 
detailed knowledge of the detector properties, which is out of the reach of a theoretical analysis. 



See also Ref. 59 for a study of gg — > hh in composite Higgs models, and Ref. 60 for an analysis 
of this process beyond the parton level in the context of the SM. 

Additionally, we studied if applying a cut on the invariant mass m^u could be useful to measure 



deviations from the SM cross section. Therefore, in Fig. 12 we show the same as Fig. 10 but after 
an invariant mass cut of mhh > 600 GeV has been applied. As can be inferred from the plot the 
composite cross section is more enhanced compared to the SM than without application of a cut, 
see also Fig. [10| On the other hand the absolute value of the cross section after cuts becomes very 
small. The plots reveal, however, another interesting feature. While for masses of the lightest top 
partner above 2 TeV the total cross section is reasonably well approximated by the cross section 



where only Higgs nonlinearities are considered, see Fig. 10, this is not the case any more after 



application of cuts. This can be inferred from Fig. 12 by comparing the full result, given by the 
points, to the black line, which is the ratio of the double Higgs production cross section considering 
only Higgs nonlinearities to the SM cross section (the full top dependence has been included in 
both cases). So we conclude that not only the heavy top partner limit in the total cross section of 
double Higgs production is a rather bad approximation unless the top partners are really heavy, 
but this approximation becomes even worse when a cut on mhh is applied. The latter, however, 
may be relevant in the experimental analyses to enhance the signal to background ratio and to 
extract information on the couplings involved in the process. 

7 Conclusions 

Models of electroweak symmetry breaking aiming at giving a rationale for the stability of the weak 
scale under radiative corrections predict an extended top sector at an energy scale typically below a 
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TeV. Carrying color and electric charges, these top partners are naively expected to give significant 
corrections to the loop-induced couplings of the Higgs boson to massless photons and gluons. We 
examined this question in the context of composite Higgs models where the Higgs boson emerges 
as a Goldstone boson from a strongly-interacting theory. We first extended the well-known SM 
Higgs low-energy theorem that gives a simple way to estimate the contribution of heavy particles 
to the Higgs couplings to photons and gluons and we then checked the accuracy of this LET 
approximation to an explicit full one-loop computation taking into account the contributions of 
all fermionic resonances. We confirmed that in composite models there is actually a quite efficient 
cancellation for the contribution of the top partners to the single Higgs production cross section, 
which deviates by no more than a few percents from the result obtained taking into account the 
Higgs nonlinearities only. For single production, the LET provides a very accurate prediction of 
the cross section. The situation is, however, totally different for double Higgs production for which 
the LET approximation is not reliable any more, and deviates from the true result by more than 
50%. The top partners also significantly reduce the enhancement of the gg —> hh cross section over 
the SM that was previously computed taking into account the strong dynamics effects only. This 
dependence on the top partner spectrum and couplings gives an indirect access to this sector that 
will complement the information gathered from direct searches as well as from electroweak and 
flavor precision data. 

The imminent discovery of the Higgs boson will put the identification of its true nature on the 
immediate agenda of high-energy physics and a careful study of the top sector can bring invaluable 
information. If the Higgs is a composite object, then tih and gg — > hh will be important channels 
that can give access to the top partners and where large deviations compared to the SM predictions 
are expected. 
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Appendices 

A Derivation of the hgg, hhgg and ^77 couplings in the SILH formalism 

In this section we derive the expressions of the couplings hgg, hhgg and /177 in the SILH formalism. 

A.l The hgg and hhgg couplings 

Our starting point is Eq. Q. We also stress that since we are working in a general basis where 
c r 7^ 0, the relation between (H) and v is non-trivial, as can be read off the W boson mass term, 




with 




(69) 
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We assume for definiteness the presence of one or more vector-like top partners, which upon inte- 
gration contribute to c g , and identify the light mass eigenstate of the heavy fermion mass matrix 
Ai with the top quark, whose mass reads 

where Cy parameterizes the correction to the SM top Yukawa coupling. The coefficients Aip in 
Eq. ^ can be related to Cy and c g by separating the contribution (to the hgg and hhgg coupling, 
respectively) of the top quark, which involves c y , from that of top partners, which is parameterized 
by c g . The results are 

1 . If d 



(H)A 1 



\(H) 2 A 2 



dlogH 
d 2 



logdetM 2 (H) 
d 



H=v 



f 



Vt 2 
o 3Cq ttV 

nip 



dlogH 



logdet M 2 (H) 



■1 



-(*) 



H=v 



f 2 mj 



(71) 
(72) 



^d(logtf) 2 

where we work at 0(1/ f 2 ) 13 Note that the 'implicit' expressions containing the determinant are 
in practice more useful than the explicit ones written in terms of Cy ^ and c g , because using the 
former avoids diagonalizing the heavy fermion mass matrix, a rather complicated task in presence 
of multiple top partners. 

In Eq. ^ we have assumed that h has a canonical kinetic term. However, in the SILH La- 
grangian the operators proportional to ch and c r correct the Higgs kinetic term as follows 



AC 



h kin 



1 

2p 



r„ +-)((H) + h)%hd»h, 



(73) 



which also contains Higgs derivative interactions. At first order in £, these can be eliminated by 
the nonlinear redefinition f8l 



h-*h- 



CH + 



h — —\ 

v 3v 2 J 



(74) 



which leaves h canonically normalized. Notice that in a nonlinear cr-model, the Higgs is canonically 
normalized at all orders, which corresponds to the relation ch = — c r /4. Performing the transfor- 



mation in Eq. (74), we arrive at the effective coupling of the Higgs to one and two gluons, Eqs. (Ill 



and (12), respectively. The invariance of these expressions under the reparameterization in Eq. ¥1 



can be easily verified by using Eqs. (71) and (72), respectively. 



A. 2 The /177 coupling 



Starting from Eq. (15), recalling the expression of the W boson mass in Eq. (69) and taking 



into account the rescaling needed to make the Higgs kinetic term canonical, see Eq. (74), it is 



straightforward to obtain Eq. ( 16 ). Similarly to Eq. ( 71 ), A\ can be related to Cy and c 7 (assuming 
that all contributions to c 7 come from the fermion sector) 



1 



(75) 



2 N ' 1 2\3\ogH v ')^.. y f 2 'm 2 2Qi 

Plugging Eq. (75) into Eq. (16), the invariance under the reparameterization in Eq. ^ is explicit. 



(H) and v expressed by Eq. ( 69 1 is higher order in £ there 



In the second and third term of each of Eqs. ( 71 1 and ( 72 1 we have used the fact that the distinction between 
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B The SILH coefficients for the Littlest Higgs and MCHM4 



In this section we give the minimal details which are needed to compute the coefficients of the SILH 
Lagrangian relevant to Higgs production via gluon fusion in the Littlest Higgs and in the minimal 
composite Higgs model with fermions in the spinorial representation. 



B.l Littlest Higgs 



The Littlest Higgs model [61] is based on the coset SU(5)/SO(5). We consider here a variation 
where only one U(l) group is gauged, as this eliminates one source of custodial breaking and thus 
relaxes the tension with EWPT suffered by the original model. In Ref. [62] it was shown that a 
scale as low as / ~ 1.2 TeV is allowed in this case. This, however, leaves an extra singlet Goldstone 
boson in the spectrum, whose effects will be ignored in the following. 14 The S field reads 



£(x) 



e 2in ^E , 



/r//(2 v / 5) H/y/2 ip 
n= Hi/y/2 -2?y/V5 H T /y/2 
\ H */V2 ri/(2V$). 
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)■ E H 
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(76) 



where H is the Higgs doublet, (pis & complex triplet and r\ is a singlet. An SU(2)\ x SU (2)2 x U (l)y 
subgroup of the global symmetry is gauged and is spontaneously broken at the scale / to the 
diagonal SU{2)l x U(l)y. In Eq. (76) we omitted the GBs that get eaten by the heavy SU(2) 
triplet of vectors. The two-derivative Lagrangian reads 

f 



C 



Tr 



D M E = ^E-t 9jWj(QjY; + EQ" T ) — ig'B^VE + TY) , (77) 

i=i,2 



with the gauged generators given by 
f a a /2 



Qi 



Q a 2 



-a a *j2 / 



y = diag (1/2, 1/2, 0,-1/2, -1/2) . (78) 



The SM fermions are taken to transform only under SU{2)\ x U(l)y- 

The coefficients ch and c r receive contributions of three different kinds. The first arises from 
the nonlinear a model structure. Using 



g 2 f 2 sm< 



\V2fJ ' 
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9192 
V9i + 92 



(79) 



we find c a H = 1/6 and c% = -4c^ = -2/3 . 

The second contribution comes from integrating out the heavy vector triplet. The procedure 
has been described in detail in Ref. [27|, and we simply apply it to the case under study, obtaining 
c v H = 1/4 and c" r = — 1 , in agreement with Ref. (To) . 

The third and last contribution arises from integrating out heavy scalars. Since we also need to 
compute eg, we write down the scalar potential up to order H 6 . Neglecting g' , the relevant terms 
are 
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'Additional sources of symmetry breaking are needed in order to give the singlet a potential. 
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where the coefficient c+ receives contributions from g\ , whereas c_ from g 2 and from the top Yukawa 
sector. In general, starting from a Lagrangian of the form 



£ $ = -$ a *n<f> a -(M' z -p 2 \H\ 2 )<!> a *$ a +( /3/$ a *ir T ey# + h.c. J + ^j-^> a *H T e^-H\H\' z + h.c. 

(80) 

(e = io 2 ) and integrating out one obtains c s H = /3 2 / 4 /(2M 4 ) and c s r = 2/3 2 / 4 /M 4 . In addition, 
there is a contribution to cq, 

c 6 = -T + ^717^ . ( 81 ) 



A \ M 2 2M 4 

where A is the Higgs quartic coupling (after the triplet has been integrated out). In the Littlest 
Higgs case we make the identifications 

M 2 = (c + + c_)/ 2 , /3 = ~(c_-c+), /3 2 = ^( c+ + c _), /3 4 = ^(c + -c_). (82) 



Therefore we find 



1 / c_ - c+ \ 2 s / c_ - c+ x 2 



*-ilf+£J ■ ^-Ifr?' ' (83) 



1 Z' c+ + c_ \_ , $ 2 



On the other hand, 

«. = -H^r=J +<* = -§■ («) 

where we have used the expression of the quartic coupling A = c+c_/(c + + c_) . Notice that in 
general the neutral component of cp gets a nonzero VEV, which is strongly constrained by EWPT. 
Small values of / ~ 1 TeV in fact require the approximate condition c+ ~ c_ to be satisfied, which 
makes the triplet VEV very smal 15 62 . We assume this condition to be realized, and therefore 
neglect effects due to the triplet VEV in our discussion. 

Concerning the top sector, in addition to the doublet qL = {^l-,^>l ) T and to the singlet £r a 
pair of SU (2)-singlet fermions Tl, Tr with electric charge Q = Y = 2/3 is introduced. The Yukawa 
Lagrangian then reads 

- £y = ^ftReijkeabXi^ja^kb + A 2 /TrTl + h.c. (85) 

k = 1, 2, 3 and a, b = 4, 5). Here x is an SU(3) triplet, x = ( &t > > ) T - The fermion mass 
matrix in the Higgs background 

* = Iff/Vl < 86) 



reads 
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implying 
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detM f M = -A? A 2 ,/ 4 sin 

2 / 



15 When c+ = c_ the potential does not contain any tadpole for 
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This allows us to write the amplitude for gg — > hh in the Littlest Higgs, in the low-energy theorem 
limit as 



Clet("S) 



s — ml 



3 \ c_ + c + 
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C_ — C4 



C_ + C i 



(89) 



We also note that in this case it is easy to diagonalize explicitly the fermion mass matrix at 0(1/ f 2 ), 
obtaining 

-C m = m t (H) t R t L + m T (H) T R T L + h.c. (90) 

where 



m t (H) 



V2 



1 + 



H 2 



1,1 
'3 + 4A T 



with = \/2 A1A2/1/AI + A| and At = \/Af + A| • In the notation of Eq. (20), we have 
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Thus — 2clt = 2/3 = const., as must be the case since the factorization in Eq. (19) holds. 



B.2 MCHM4 

Similarly to what we did for MCHM5, we can apply directly Eq. ^ to derive the hgg and hhgg 
effective couplings at all orders in £. For more details about the model we refer the reader to 
Ref. (9], the notation of which we adopt here. The composite fermions are embedded into SO (5) 
spinors 4 ly / 6 as 



ipi 



QLl 
-iT L 
-iB L/ 



Qr 
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(94) 



where q L \ = (t L i,b L i) T and Q R = (T R ,B R ) T are SU(2) L doublets, and T L , B L ,t R1 ,b R i are 
singlets. In addition, an elementary doublet qii = (*L2j&L2) T and singlet t R 2 are present. The 
fermion Lagrangian reads 



C f = yi> L T M $ M ^R + fX q q L ,Q R + fX t T L t R2 + h.c. , 



(95) 



where 3> = /(0, 0, 0, sin(H/ f), cos(H/f)) , and T M are the Gamma matrices of SO(5) . The mass 
matrix reads 



tR1 

~Cm = (t L i t L2 T L ) M \t R2 I !- h.c. 

\ T R, 

from which we find 
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detM^M = X 2 q X 2 t y 2 fsm 2 (H/f) and At = * y/l - £, A 2 = ~ . 



(97) 
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Finally recalling the expression of the Higgs trilinear coupling in MCHM4, £^3 = —(mf l /2v)h' i y/l — £ , 
we find the amplitude for Higgs pair production via gluon fusion in the low-energy theorem ap- 
proximation 

C%g HM \s) = T M^(l-0-l. (98) 
s — mi 



C The x 2 test for electroweak precision observables 

We discuss here in detail the \ 2 test used to constrain the parameters of the MCHM5 described 



in section 4.1 The best experimental determination of e%, 63 and e& still comes from the precision 
measurements at the Z pole mass at LEP 
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Here p is the correlation matrix between the ej obtained from the App. E of Ref. 39 , marginalising 
over the three parameters Mz, as(Mz) and Aa\' d (Mz)- The status of electroweak precision 
observables did not change since then, except for the mass of the W boson. The latter was recently 
updated based on Tevatron results 1 63 , 64 1 , and the new world average is now [65] 



m w = 80.385 ± 0.015 GeV. (100) 
The parameter €2 is the only one depending on the mass of the W boson, through the term 
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where sq = 0.23098, and the measurement of Ar w is related to mw through 
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Here a(0) is the fine-structure constant and Gf the Fermi constant, both known to high accuracy. 
Furthermore, Aa accounts for the running of the electroweak coupling between the low energy limit 
and the Z-pole mass. The uncertainty associated to it is important, but the shift in Ar w induced 
by the new value of mw is independent of Aa. The change of mw and subsequently €2 between 
the LEP data of 2006 and the present value is then 
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The experimental values for the used in this paper are therefore 
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where the changes in the second row and column of the correlation matrix p come only from the 
reduction of the uncertainty on £2; the covariance matrix between the ej remains unchanged. 
On the theoretical side, the e% are predicted to take the values j24j 



ef l) = [+5.60 - 0.86 log (m h /m z )] • lO" 3 + Aef + Ae{™ ns 

4* fc ) = [-7.09 + 0.16 log (m h /mz)]-W- 3 

e^ h) = [+5.25 + 0.54 log (m h /m z )] ■ 10~ 3 + Ae^ R + Ae^ v 

e (th) = _ 6 . 43 . io-3 + Ae ^™ ns 

where the first numbers are the Standard Model corrections and the remaining contributions are 



given in Section 4.1 for the MCHM5. For the computation, we used a Higgs mass of = 125 GeV. 
The x 2 test is then defined as 

X 2 & d> L , *r, R, m p ) = £ {4 h) ~ zf^) (ef > - P) ) , (106) 

id 

where C _1 is the inverse of the covariance matrix 

C ij = Ae < T p) Pi j Ae ( f xp \ (107) 



As indicated in Eq. ( 106 ), the x 2 depends on five parameters. However, £ carries a different meaning 
than the other parameters, since it provides a measure of the fine-tuning of the model. The absolute 
minimum of the \ 2 is in particular obtained for a very small value of £, which is highly unnatural. 
We compute therefore the minimum of the x 2 f° r a fixed value of £, and require 

X 2 (<pL,<f>R,R,m p )\ ( - Xminlt ~ 13 ' 28 ' (108) 

The value on the right-hand side corresponds to a confidence level of 99% with four degrees of 
freedom tpR, R and m p ). For the values quoted in section 4.1, we have 



Xmm|^ = o25 0.98, Xmin|< = Q .1 0.85 . (109) 

Note that the minimum of the x 2 in the MCHM5 is significantly lower than in the Standard Model, 
X%m = 4-84, which is expected due to the larger number of fitting parameters. 

D Partial decay widths of heavy fermions in MCHM5 

In this section we collect the formulae for the partial decay widths of fermionic resonances in 
MCHM5. We start by defining the relevant couplings. We denote by Ul,r the transformations 
that diagonalize the mass matrix in the top sector, 

M ->• Ul M Ur = diag (mi, m 2 ,7713,7714) , (110) 



where A4 is the mass matrix in the basis where the rotations in Eq. (45) have already been per- 
formed 16 On the other hand, G h jj is the Yukawa coupling matrix after the rotations in Eq. (45), 
and R are the matrices containing the couplings of the fermions to the Z boson, 

G Z L =diag (-- . - - -4 . - - - 4+ . -4+ j . (Ill) 
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16 The mass eigenstates fiL,R (i = 1, . . . , 4) are ordered by decreasing mass, so Ja is identified with the top, and fs 
with the lightest top partner. 
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for left-handed and right-handed fields, respectively (the ordering of the fields is understood to be 
that of Eq. (44)). Notice that the rotations in Eq. (45) leave G£ r invariant, because they only mix 



states with the same EW quantum numbers. 



D.l Charge 2/3 states 

We have for the lightest top partner ip the partial decay width into a Z boson and a top quark 
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where gz = g/ cos6 w . From the Yukawa Lagrangian Eq. (57) we can extract the leading- 
couplings of T, T, X 2 /^ . They are given by 
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where Xwl = (9/V^)(Ul) 13 • The leading order couplings read 
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For the decay tp — > ht we find 
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D.2 B and X 5 / 3 

For x = A 5 / 3 , B we find for the decay width into a W boson and a top quark 
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The leading order couplings are 
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The leading order expressions of the heavy fermion masses are reported in Eq. ( 46 ) 



The formulae for the partial decay widths contained in this Appendix reduce to those given in 



Ref. |66| when the approximations in Eqs. (116, 117 119 123 124 129) are made. 



E Analytical results for the gg — > hh cross section in MCHM5 

We present here the analytical result for the partonic gluon fusion cross section into two Higgs 
bosons, dgg^hh, in MCHM5. 

E.l Notation 

The four-momenta of the gluons are denoted by p\ and p2, and the four- momenta of the Higgs 
bosons by p 3 and p^. All momenta are taken incoming. The Mandelstam variables s,t,u are given 
by 



8 = (Pi +P2) 

The scalar integrals are defined as 

ri I 2 2 2\ 

Ci j (m 1 ,m 2 ,m 3 ) = 
d 4 q 



t = (Pi +P3) 



U = (p2 + P3) 



ivr 2 ( g 2 _ m 2) (( ? + p .f _ m 2^ q + p . + Vi) 2 _ m 2) 

D ijk (m 2 ,rn2,ml,ml) = 



(130) 



d 4 q 



1 



\2 2\ 



ivr (g 2 -mf) ((g + ft) - m 2 )((q + Pl + Pj ) - m 2 )((q + Pi + Pj + p k ) - „ H) 
The analytic expressions can be found in Refs. [671 . They have been evaluated numerically in our 



code with the help of Loop Tools 68 
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E.2 Tensor basis and projectors 



The tensor basis, which has been given in Ref. 13 , reads 



V I 1 

P1P 2 

(Pi -P2) 

J2 „v „M 



P\P\P 2 2(P3-P2)PiP3 ^{PZ-Pl)P Y iP2 , 2 PsP3 



and 



P%{P\-P2) Pl(Pl-P2 

with 
Ai ■ A 2 = 



pl {pi -P2) 



+ 



Pt 



(131) 
(132) 



2 _ (glVP3) (P2 -P3) 2 



Pt - 
and 



(Pi -P2) 
Ai ■ A! = A 2 ■ A 2 = 2 



(133) 



E.3 Triangle form factor 

The triangle form factor can be cast into the form 

Fjf\(m) = 2 [2m + (4m 3 — s m) C\ 2 {jr? , m 2 , m 2 )] 
and can be found in Ref. 



(134) 
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In the limit of large quark mass m S> \/S ~ rrih the triangle form 
factor simplifies to F/\ = 21/ (3m). This is equivalent to applying the low-energy theorem. The 
corresponding amplitude is given by 



Aa = Y-Al u e*e b J ab ^2 ( <■ _ X 2 9hhh 9hq m +2g hhMi ) F A (rm 
47r i= i V s m h / 



(135) 



The couplings ghq^i and ghhcjiqi are the diagonal elements obtained from the Higgs coupling matrices 
y Ghff in Eq. (49) and y/(2 /) G hh ^j in Eq. (50 ), respectively, after rotation to the mass eigenstate 



basis. The tripleHiggs coupling ghhh is given in the MCHM5 by 

2£ 



ghhh 



3ml 1 



(136) 



In the SM limit, in Eq. ( 135 ) there is no sum over heavy top partners, and we are only left with the 



top quark contribution with the Higgs coupling to the tops given by g^tt = mt/v, and gthtt = 0. 



The triple Higgs coupling in the SM limit can be obtained from Eq. (136) by setting £ = 0. 



E.4 Box form factors 

In the box diagrams we can have spin S z = and 2 gluon gluon couplings. The matrix elements can 
therefore be written in terms of two gauge invariant form factors. Furthermore, we have divided 
the form factors for the boxes into the parts which do not involve a 75 and the parts which are 
proportional to two 75 matrices. Couplings with a 75 arise only for Higgs couplings to two different 
fermions (but of same flavour). The diagrams including only one 75 vanish because of the sign 
flip of the coupling when the direction of the fermion line changes. The form factors have been 
calculated with FeynCalc 



69 



and checked against Ref. |70|. For the limit mi — > m 2 they agree 
with Ref. [14]. The form factors are UV- finite as the coefficients in front of the UV-divergent one- 
and two-point functions Aq and Bq are anti-symmetric in rm and rrij and vanish upon summation 
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over i, j. 

We introduce the following abbreviations 



Cl2 = 


ri { 2 2 2\ 


C13 


r< 1 2 2 2\ 


C14 = 


Ci4(m±, m 2 , m 2 ) 


C23 


= C 23 {m\,m\,m 2 2 ) 


C*24 = 


ri 1 2 2 2\ 
C24(m ll m 1 ,m 2 ) 


C34 


ri 1 2 2 2\ 

= G u (m 1 ,m 1 ,m 2 ) 


D 123 = 


n / 2 2 2 2\ 


-D132 


7-1 / 2 2 2 2 
= £>i32(m,i, m 1 ,m 2 ,m 2 


-D213 = 









(137) 

The box form factors Fn and G*n associated with spin and spin 2, respectively, are then given by 
(contributions which cancel by summing up all contributions are omitted) 



Fn(mi,rrij 
2 r 



2s + \m i s C\2 + s((m, + nij){2m i (m; + rrij 



mis) - m 2 (i + u))(D 123 + D 132 + D 



+ (m 2 , - (rrn + mj) 2 ) [(t - m 2 h )(C 13 + C24) + (fi - m^Caa + C U ) 
- (iu -m\ + s(m 2 - m 2 ))D r32 ] 



213 J 



:i38) 



Gn{mi,mj) 
1 



i-u — 



(P + - (4m^ + Amimj)(t + u) + 4(m 3 - - m i )(m i + m^) 3 + 2m\)sCi 2 



+ (m\ + 1 2 - 2t{mi + m i ) 2 )((t - m^)(Ci 3 + C 24 ) - stD 213 ) 
2u{nii + mj) 2 )((u - ml)(C 23 + C u ) - suD 123 ) 
(P + u 2 - 2m^)(t + u - 2(m» + m,) 2 )^ 

(t + u - 2(m; + m J ) 2 )((tu - m£)(ro? + m 2 ) + s{m 2 - m 2 ) 2 ){D 123 + D132 + #213) 



'72 



and 



F n ^(mi,mj) = -F u (mi,-mj), Gn^im, m>j) = -Gnim, -mj). (139) 

Here i*h 5 and Gn,5 denote the spin and 2 box form factors which are proportional to the Higgs 
couplings to quarks, containing a 75 matrix. In the large quark mass limit for rrn = Wj the 

form factors reduce to Fq = —2s/(3mf) and Gn = Or^l The spin and spin 2 box amplitudes read 



44 44 

i=l j'=l i=l j'=l 



and 



44 44 

i=l j=l i=l j=l 



For rrii — rrij the couplings in front of the form factors F\j,s and G\j,s vanish, so that in this case these form 
factors are not needed. 
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The couplings g h q iqj and ghq iqj ,5 are given by 



Qhqiqj 
9hq~iqj,5 



G hff,ij + G h ffj % 
G hff,ij ~ G hff,ji j ■ 



(142) 

(143) 



where Gh / fij denotes the (ith, jth) matrix element of the coupling matrix of Eq. ( 49 ) in the mass 
eigenstate basis. Note that ghq iqj £ is antisymmetric in i and j and hence changes sign, if incoming 
and outgoing fermions in the vertex are interchanged. 
The complete amplitude of the process is given by 



A(gg -> hh) = A A + A ,n + A 2 ,u 



(144) 



E.5 Expansion of the form factors 

In the case where top partners are neglected, we can perform explicitly the expansion of the form 



factors in 1/m 2 , i.e. for small external momenta (see Ref. 71 ), and go beyond the leading order, 
the latter corresponding to the LET result. We find 



^A 
Fa 
Go 

The partonic cross section then reads 



+ 



7 S 



m t \3 180 mf 
air \ 3 30 m\ 



s 11 / m 



Vgg^hh 



1024(2vr) 3 s 2 J { 



mf 90 



dt 



tu 



(145) 
(146) 
(147) 



\C A F A + C n F n \ 2 + \C n G n \ 2 



with t± given by Eq. ( 68 ) and 



\C A F A + CuFuV + \CuGuV = ^(ca - cti) 



, 2 = £4. 
v 4 9' 



1 + 



1 



CA 



3ml 



1 - 2£ 



4e, 



m^CA - cn 

Cn = {7T^ 



7 - 7 2 



60 

2 



10 



where we used 



s — mt 



-Qhhh 9htt + 2 5h/iit 



mt 

~YC A , 



-cn 



(149) 
(150) 

(151) 



with the couplings given by g Mt = (m t /v)(l - 2£)/Vl - £, g hhU = -2m t £/v 2 and g hhh as defined 
in Eq. (136). The leading term in 1/mf corresponds to the LET result, see Eqs. (28) and (66). 
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